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ABSTRACT 

We examine scaling relations of dispersion-supported galaxies over more than eight orders of magnitude in 
luminosity by transforming standard fundamental plane parameters into a space of mass, radius, and luminosity. 
The radius variable ri/2 is the de-projected (3-D) half-light radius, the mass variable Afi/2 is the total gravi- 
tating mass within this radius, and L1/2 is half the luminosity. We find that from ultra-faint dwarf spheroidals 
to giant cluster spheroids, dispersion-supported galaxies scatter about a one-dimensional "fundamental curve" 
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through this MRL space. The mass-radius-luminosity relation transitions from AI1/2 

faintest dwarf spheroidal galaxies to Afi/2 ~ rl'l^^ ^ ^I'A ^'^^ '■^^ most luminous galaxy cluster spheroids. 
The weakness of the M1/2 ~ ^1/2 slope on the faint end may imply that potential well depth limits galaxy 
formation in small galaxies, while the stronger dependence on L1/2 on the bright end suggests that baryonic 
physics limits galaxy formation in massive galaxies. The mass-radius projection of this curve can be compared 
to median dark matter halo mass profiles of ACDM halos in order to construct a virial mass-luminosity re- 
lationship (^fvii ^ L) for galaxies that spans seven orders of magnitude in Afvii- Independent of any global 
abundance or clustering information, we find that (spheroidal) galaxy formation needs to be most efficient in 
halos of A/vir ^ 10^^ Mq and to become inefficient above and below this scale. Moreover, this profile match- 
ing technique for deriving the A/vir — L is most accurate at the high and low luminosity extremes (where dark 
matter fractions are highest) and is therefore quite complementary to statistical approaches that rely on having 
a well-sampled luminosity function. We also consider the significance and utility of the scatter about this re- 
lation, and find that in the dSph regime observational errors are almost at the point where we can explore the 
intrinsic scatter in the luminosity-virial mass relation. Finally, we note that purely stellar systems like Globular 
Clusters and Ultra Compact Dwarfs do not follow the fundamental curve relation. This allows them to easily 
be distinguished from dark-matter dominated dSph galaxies in MRL space. 

Subject headings: galaxies: dwarf ; Galaxies: Local Group ; galaxies: elliptical and lenticular, cD ; galaxies: 
fundamental parameters; cosmology: dark matter 
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1. INTRODUCTION 

Galaxy observables such as size, luminosity, and velocity 
dispersion are known to follow scaling relations. The study 
of these relations provides a window into the processes that 
regulate galaxy formation. The ACDM dark matter halos that 
host these galaxies are also predicted to follow structural scal- 
ing relations, including relations between their central densi- 
ties and total virial masses. In this paper, we seek to link 
galaxy observables to dark matter halo properties by studying 
galaxy dynamical masses (M1/2) within their 3-D half-light 
radii {ri/2) ^s a function of galaxy luminosity (L1/2 == L/2). 
This coordinate space of intrinsic parameters (MRL Space) is 
obtained via a simple transformation of the standard observed 
parameters of fundamental pl ane space. Our approach is mo- 
tivated by the work of lWolf etalJ (2010). who showed that 
the dynamical mass of a spheroidal galaxy within ri/2 can be 
determined accurately from observed sizes and velocity dis- 
persions without knowledge of the stellar velocity dispersion 
anisotropy. This fact enables manifestly apparent physical in- 
terpretations of MRL space and, in principle, a method to con- 
nect central galaxy densities to global dark matter halo prop- 
erties. 



' Center for Cosmology, Department of Physics and Astronomy, 4129 
Frederick Reines Hall, University of California, Irvine, CA, 92697, USA; 
etolleru@uci.edu, bullock@uci.edu, wolfj@uci.edu 

^ Miller Fellow, Department of Astronomy, 601 Campbell Hall, Univer- 
sity of California, Berkeley, CA 94720; graves@astro.berkeley.edu 



It is well established that when placed in a parameter 
space of observed velocity dispersion (a), 2-D effec- 
tive radius (Re), and surface brightness (le), bright 
(> L*) early-type galaxies lie appr oximately within a 
two-dimensional "fundamental plane" (Diorgovski & Davisl 
1987; Dressier 1987; Faberetal. 1987). Other work (e.g . 
Nietoetal. 1990; Bender etal. 1992; Burstein et al. 1997t 
Prada & Burkert 2002; Zaritsk v et al 2006b: Shankar et all 
booet iWoo et al. 2008; Forb es et alJ l2008h has expanded 
upon or considered similar such relations sometimes in- 
cluding galaxies that have significant rotationally-supported 
components. These scaling relations provide a wealth of 
opportunities to examine wh at physical proce sses gen- 
erate them (e^. Dantas et al J 120001: jDekel & W oo 20031 

all 120081; iHopkins et aD 



Robertson et al.' 2006; ZaritskiLet 



2008 



Kormendv et al. 2009|; |Boyill & Ricotti 2009 



Graves et alj|2009b) . and hence further constrain scenarios of 



galaxy formation. 
Zarisky 



and 



collaborators 
(IZaritsky. Gonzalez. & Zabludofil l2006allbh explored a 
unified description of the fundamental plane parameters 
for all spheroids that are embedded within their own dark 
matter halos. They found that dwarf spheroidal galaxies 
(dSphs), dwarf elliptical galaxies (dE), normal elliptical 
galaxies (E), and the extended stellar spheroidal components 
of galaxy clusters (cluster spheroids, CSphs) could be 
characterized by a 2-D fundamental manifold in (a, Re, 
le) space (although curved relations had been seen noted 



in sub-spaces, e.g. lGrahaml|2005h . The CSph of a galaxy 
cluster halo is the sum of the brightest cluster galaxy (BCG) 
and the extended intra-cluster stars (ICS). Empirically, 
the inclusion of CSphs is motivated by the fact that they 
demonstrate a relationship between R p and If, tha t is similar 
to elliptical galaxies in many respects (iGonzalezet al. 2005) . 
From a theoretical/cosmological perspective, a CSph is the 
most natural single stellar system to associate with the host 
dark matter halo of a cluster (while the cluster galaxies 
themselves are more readily associated with subhalos, e.g. 
IConroy et al. 2007; Purcell et al. 2007). Typically, the CSph 
(L > lO ^^Lp,) contains many more stars than the BCG 
by itself JGonzalez et alj l2005h and, importantly for our 
purposes, it extends to a fair fraction of the cluster virial 
radius, and thus (in principle) allows a more global probe 
of the cluster potential. We will include CSphs as the 
cluster-halo counterparts to normal spheroidal galaxies in our 
work. Hence, our definition of "galaxy" here is the central 
luminous component of a distinct dark matter halo (although 
it may be a subhalo of a larger halo, as is the case for dSph 
satellites or non-BCG galaxies in clusters). 

On the opposite luminosity/size extreme from giant CSphs 
are the ultra-faint dw a rf dSphs (e.g. [Wil l man et al. 2005|; 
Belokurovetin 120071: lWalshetaljl2^ iBelokurov et al] 



20091) . Discovered by searches within the Sloan Digit al Sky 



Survey dataset (SDSS. lAdelman-McCarthv et al.l2008l) . these 
systems can have luminosities smaller than ~ IOOOLq and 
have been shown to be the mos t dark matter dominated sys- 
tems known (e.g. Martin et al.' 2007; Simon & Geha 2007; 
Fenarrubia et al.'2008b; Geha et al. 2009; Simon et al. 2010) . 
The ultra-faint dSphs provide means to study galaxy for- 
mation within the smallest dark matter halos that host stars 
JStrigari et al.l l2008b). By including them in our analysis, we 
will extend galaxy scaling relation studies to span more than 
eight orders of magnitude in luminosity. 

By exploring galaxy properties over a very wide range in 
luminosity, we are able to address one of the broader ques- 
tions in astrophysics: how and why galaxy formation effi- 
ciency varies as a function of dark matter halo mass. Re- 
markably, observed galaxy luminosity functions and two- 
point clustering statistics can be explained fairly well under 
the assumption that L (or stellar mass) maps to dark matter 
halo mass Myjr in a monotonic way (e.g.lKraytsov et al.l2004l : 
IConrov & Wechsleiil2008l: iMoster et alj|2009l) . such that halo 
viral mass-to-light ratio M^i^/L is minimized near L ^ L^, 
and rises steeply at larger and smaller L. An understanding of 
this behavior - how and why it happens - is hampered at the 
smallest and largest mass scales because luminosity functions 
become less complete and less well-sampled in the extremes. 
One of the goals of this work is to use galaxy MRL relations 
to inform the Mvir — L mapping in a way that is independent 
of large-scale abundance and clustering studies. 

The paper is organized as follows: in i|2] we describe the 
data set used for this study and the relevant sources. In ij3] 
we consider the scaling relations of our data, introducing a 
new space ("MRL" Space) that is designed to provide a bridge 
between the scaling relations of galaxies and the scaling re- 
lations of dark matter halos over the full dynamic range of 
known galaxies (^ 10 orders of magnitude in L and Afvii). 
In 21 we introduce a one-dimensional curve that the galaxies 
follow in this space, and apply this curve t o can onical ACDM 
halos to map halos onto their galaxi es. In i l4.3l we address the 
scatter in the fundamental curve, in il5.2l we address the errors 
and scatter in the halo mapping, and in [j6] we conclude. 



Throug hout this paper we assu me a ACDM cosmology with 
WMAP7 (iKomatsuet all 120101) parameters of h = 0.704, 
Qm = -272, r^A = 1 - ^Af, CTs = 0.809, and n, = 0.963. 
Further, we use the symbol log to represent base- 10 loga- 
rithms. 

2. DATA 

The data sources for this study are varied by necessity due 
to the wide dynamic range covered. Table [T] gives the rele- 
vant parameters for the objects in this study and the sources 
for each. Starting with the least luminous objects that are em- 
bedded within dark matter halos (L < IO^Lq), our dwarf 
spher oidal (dSph ) data set is taken from the summary ta- 
ble of I Wolf et al.l (120101) and draws from various sources for 
photometric properties and resolved star kinematic measure- 
ments for Milky Way dSph galaxies. Moving up in brightness 
(L ~ IO^^^Lq) our "dwarf elliptical" (dE) sample is taken 
from the Virgo Cluster dE study of Geha et al. (2003). Note 
that while dEs are not as clearly dark matter dominated as 
dSph galaxies within ri/2 (see below), they are believed to be 
embedded in their o wn dark rnatter h alos based on extended 
kinematic samples ( Geha et al.ll2010l) . 

Data for normal elliptical galaxies (E) are from lGraves et al.l 
( I2009ai L ~ lO^^i©) and are discussed in more detail toward 
the end of this section. The brightest (L ~ lO^^LfT)) cluster 
sphero id (CSph) data are from the imaging o f iGonzalez et al] 
( 120051) and spectra of IZaritskv et all (I2006H) . These data are 
also described in more detail below. 

We also examine two comparison populations as exam- 
ples of systems that are not embedded within dark matter 
halos: Milky Way globular clusters (GCs, L ~ 10^) and 
ultra-compact dwarfs (UCDs, L ~ 10^). For GC photometr 
we use the 2003 revision of the Ha ms catalo g (Harris 199' 
and take velocity dispersions from Prvor & Meylanj fll993b . 
For UCDs we use data fromiMieske et al. (2008). Note that 
while the status of UC Ds as large examples of purely stellar 
systems is debated (e.g.lEvstigneeva et al .1120071; iGoerdt et ajj 
2008; Baumgardt & Mieske 2008; Dabringhausen et al. 200i^ 
iTavloret al.ii2010 , and references therein) we find that their 
scaling relations are more in line with GCs than similarly lu- 
minous dSphs and therefore treat them as lacking dark matter 
halos below. 

The CSph data set stands out compared to the other data 
sets in two distinct ways. First, while all other d ata sources are 
in the V band, the CSph data (summarized by IZaritskv et al.l 
'2006b) use Cousins I-band luminosities. We convert these 
data t o V-ban d using averaged colors of E galaxies from 
Fuku gita et al.l (ii995). While this does not account for the 
possibility of a systematic error in R^ for the CSph data points 
due to a different choice of band, this effect is likely to be 
small given the large dynamic range in this data set. Further- 
more, La Barbera et al. ( 2008) find that the fundamental plane 
for early-type galaxies is nearly independent of band from op- 
tical to the K-band. Given the similarity of the stellar popu- 
lations for those galaxies and the CSph, it is therefore likely 
that the band mismatch is not a significant effect. 

The second way that th e CSph data set differ s is that the 
velocity dispersions from (Zaritsk y et al.ll20"06bh are derived 
from galaxies in the cluster, rather than the CSph (mostly 
ICS) light itself. This is of course not ideal, but the measure- 
ment of ICS velocity dispersions is very difficult with current 
spectroscopic capabilit ies. While this ha s been accomplished 
both in integrated light ("Kelson et al.'2002) and planetary neb- 
ula kinematics (e.g. ,Arnaboldi et al.i i2004,) for a few clusters. 
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TABLE 1 
Key properties of dispersion-supported stellar systems. 



Name '' 


log(a) b 


log(iJe) '■■ 


log(ri/2)'i 


log{Li/2) = 


log(Mi/2)f 


log(AfDM) g 


Object Type '' 


^obj' 


Source J 


A0122 


2.83 


2.03 


2.15 


11.3 


13.7 


13.7 


CSph 




1 


E Bin 1 


1.92 


0.04 


0.16 


9.22 


9.85 


9.5 


E 


36 


2 


VCC452 


1.38 


-0.15 


-0.02 


8.04 


8.57 


8.37 


dE 




3 


Draco 


1.0 


-0.66 


-0.53 


5.03 


7.32 


7.32 


dSph 




4 


47Tuc 


1.31 


-2.66 


-2.54 


5.20 


5.92 


5.92 


GC 




5,6 


F-19 


1.36 


-1.05 


-0.92 


7.00 


7.64 


7.64 


UCD 




7 



Note. — This table will be published in its entirety in the electronic edition of the Journal. A portion is shown here for guidance regarding form 
and content. 
'^ Name of the object. 
'' Log of Velocity dispersion in km/s. 
■-' Log of 2-D Half-light/effective radius in kpc. 
■^ Log of 3-D (deprojected) Half-light radius in kpc. (see ij3) 
■^ Log of V-band half-luminosity in Lq, i.e. Ly /2. 
' Log of Half-light mass from Equation[T]in Mq . 
^ Log of CoiTected/dark half-light mass in Mq as described in 31] 

'' Object Type: CSph = Cluster Spheroid as described in Zaritskv e t alj 42006hl) . E=bright Elliptical galaxy, dE=dwarf Elliptical, dSph=local group 
dwaif spheroidal, GC=Galactic globular cluster, UCD=Ultra Compact Dwarf. 
' Number of objects per bin for E data set. 

J Source code: LZarit skv et al. 1 2006b), 2 1Graves et al.U20Q9bl) . 3 1Gehaet al.lf2003l) . 4 1Wolf et aUOOTOl) . 5 1Harrisin99a . 6 1Prvor & Mevlarj n993h . 
7 1Mieske et alj (2003) . 



there is not yet a large, homogenous sample. This is required 
for generality and to compare to our other samples, and hence 
we are forced to use galaxy dispersions until large direct mea- 
surement samples become available. In principle this could 
impose a bias in our mass estimator (described below) be- 
cause the ICS and cluster galaxies follow different distribu- 
tion functions. We e xplo re in more detail how this bias might 
affect our results in il4.3l 

The normal elliptical galaxy data comprise a sample of 
^16,00 galaxies selecte d from the Sloan Digit al Sky Survey 
(SDS S,|Yor k et al.ll2000[) Main Galaxy Samp le (IStrauss et alJ 
l2002h . as described in iGraves et alJ (12009 ah . Galaxies are 
selected to be passively evolving quiescent objects with no 
emission lines in their spectra. The individual galaxies are 
sorted into bins in the 3-D Fundamental Plane parameter 
space defined by a, R^, and 1^. Values reported here are the 
median values for each bin of galaxies. 

Before continuing, we summarize our galaxy terminology 
and the symbol codes we use when presenting each galaxy 
type. The CSph population of IZaritskv et al.l ( 1200 6b') is pre- 
sented as or ange squares. The "E " or "bright E" terminology 
refers to the iGraves et alJ (l2009al) data set and is represented 
as red circles of varying size such that the size of the data point 
is proportional to the number of galaxies in each bin. The 
dE or "dwarf elliptical" label refers to the iGeha et aP (120031) 
data set and is presented as yellow diamonds of uniform size. 
The Milky Way dSph satellites here are represented by ma- 
genta triangles. In some cases, a distinction will be drawn 
between the "SDSS dSphs" and the "classical dSphs," refer- 
ring to those discovered by SDSS and those known before. 
The SDSS dwarfs are almost exclusively fainter, and include 
the "ultra-faint dSphs." Finally, the GC and UCD populations 
are represented by the green and blue star-symbols and pen- 
tagons, respectively. 

3. MRL SPACE 

We now examine the data set described in the previous sec- 
tion in the context of the scaling relations of the observables. 
We emphasize the use of the MRL space described below to 
understand this data set. 

First, we provide a sample projection of the data set de- 



scribed in the previous section (Table [T]i- Figure [T] plots this 
data set in the 2-D space of luminosity (L) and stellar ve- 
locity dispersion dispersion {a) — the Faber-Jackson relation 
(F aber & Jackson 197 6). We also show best-fit power laws 
(L = a^) for each of our classes of objects. We compute 
slopes by fitting a linear relation in log space with log£ 
(log (j) as the parametric variable. For the CSphs, Es, dEs, 
dSphs, UCDs, and GCs, this results in slopes of 7 — 1.5 (0.5), 
2.6 (1.8), 6.0 (1.1), 11.1 (6.3), 2.4 (1.2), and 3.4 (1.5), respec- 
tively. 




6 8 

log(L/Lo) 

Fig. 1 . — Faber-Jackson jFaber & Jacksonll973) relation relating luminos- 
ity to stellar velocity dispersion for the data set. The lines show best-fit power 
laws for each of the sub-populations of objects (see text for slopes). Orange 
squai'es are cluster spheroids (CSph), red circles ai'e elliptical galaxies (E), 
yellow diamonds are dwarf ellipticals (dE), magenta triangles are Milky Way 
dwarf spheroidals (dSph), green stars are Galactic globular clusters (GC), and 
blue pentagons ai'e ultra-compact dwaii's (UCD). 



In this plane the slopes increase towards larger luminosi- 
ties, suggesting a definite scaling relation (the original Faber- 
Jackson relation). We note, however, that the dSphs, UCDs, 
and GCs are mixed together in this projection, a clear draw- 
back from interpreting these objects in this space. Further, 
there is structure to the E sample not fully aligned with this 
2-D p arameter space. The structure here is the f undamental 
plane dPiorgovski & Davis|[T987t iDressleri [19871 iFaber et"an 
Il987h for E galaxies, distinguished from the Faber- Jackson 
relation by being a 3-D parameter space with the inclusion of 
the effective radius (i?e, the radius enclosing half the total lu- 
minosity) and use of mean surface brightness /g = L/{2ttR^) 
in place of the luminosity. In Appendix iB] we show this data 
set in the fundame ntal plane space (and the related n space of 
iBender et alj 19921) for reference and comparison, but here we 
emphasize the use of a different parameter space, described 
below. 

While the fundamental plane is a valuable parameter space 
of observables, the connection to this space from typical dark 
matter scaling relations is non-trivial. In order to facilitate 
manifestly apparent theoretical interpretations, we introduce 
a set of physical variables - a mass, a size, and a luminosity 
- that are derived from the same observables. Hence, we call 
this space "MRL Space" for the three variables: 



20101 §2.1). Similar estimato rs (e. g. ISpitzedl 1 9691: llUingworthl 
19761: ICappellari et"ani2006h have the same form (by dimen- 



1. The half-light mass M1/2 = 
dynamical mass within ri/2- 



M{< ri/2) - the total 



2. The 3-D half-light radius ri/2, the radius enclosing the 
half-luminosity L1/2 ■ 

3. The half-luminosity, Li/2^ half of the total luminosity 
emitted from the galaxy {not necessarily the same as 
half the observed luminosity). 

We note that the luminosity variable here is defined in terms 
of the total luminous material in the galaxy, ignoring any at- 
tenuation that may occur as light propagates out of the galaxy. 
Below we describe the transformation of observables used to 
closely approximate this space for the data set here. 

A major motivation for the choice of these coordinates is 
the explicit use of the mass within the 3-D half-light radius 
as the mass variable, M1/2 = Af( ri/^). The adopti on of 
this mass in particular is motivated bv lWolf et al.l ( 120101) . who 
showed that while dynamical masses with r ^ ri/2 and 
r 3> ri/2 are largely unconstrained from 1-D velocity dis- 
persion data (due to weak constraints on the stellar velocity 
dispersion anisotropy), M1/2 can be determined simply and 
accurately for spherical systems without knowledge of the 
anisotropy: 

Ml/2 = 3 0-1^2^1/2. (1) 

iWolf et aTl (120 101) showed that as long as the stellar velocity 
dispersion profile is fairly flat with radius, this mass estimator 
for Ml/2 is accurate for a wide range of light profiles, includ- 
ing the types of profiles used to fit all of the types of objects 
shown in Table [T] Hence, for stellar systems with negligi- 
ble rotational support, this formula provides a good estimate 
for the total dynamical mass within ri/2 (assuming spherical 
symmetry). 

Note that Equation [T] was not derived using the virial the- 
orem, but rather follows from the Jeans Equation. The virial 
theorem provides only an integral constraint on the total mass 
traced by a stellar system and therefore cannot be used to infer 
precise masses (see Merri tt, 1987, Appendix A and .Wolf et alj 



sional analysis), but for most of these the coefficient is cali- 
brated by examining high-quality data and assuming that mass 
follows light. These calibrations are less useful for a wide 
variety of spheroidal galaxies because there is no reason to 
expect that all spheroidal galaxy are homologous. Equation 
[U is derived analytically rather than empirically, and shows 
that there is a particular radius at which the mass is unbi- 
ased at any scale (w ri/2)- Estimators that do not use this 
radius must have different virial coefficients as a function of 
scale. Further, Equation [T] assumes neither mass-follow-light 
nor isotropy, and hence is suited to the range of objects with 
various dark matter fractions that we consider here. 

Further, we note that the approximation ri/2 = 4i?e/3 is 
accurate f or the l ight p rofile s of relevance in this pa per. As 
shown in Ciotti (1991) and Lima Neto et all (11999). depro- 
jected spherical Sersic (Sersic 19631) profiles for a range of 
Sersic indicies are within a few percent of this relation , and 
the sa me is demo nstrated for Plummer (lPlurnmerlll911D and 
King dKind 1 19621) profiles in ISpitzeii (119871) and fWolf et alj 
(20101). The objects presented here are well fit by at least one 
of these profiles, motivating the use of the approximation. We 
note here that these deprojections must assume spherical sym- 
metry, like the M1/2 estimator described above. 

With these estimators chosen, the MRL space as derived 
from the observables consists of: 



1. Afi/2 = 3G' 



-1 ? 

cr ri/2 . 



'^1/2 



4i?e/3. 



3. Li/2 - L/2 - huRl 

Here R^ is the 2-D (projected) half-light radius, G is the grav- 
itational constant, a is the stellar velocity dispersion of the 
galaxy, and I^ is the mean surface brightness within R^. We 
note that the observables here are the same as those used for 
the fundamental plane and thus this space can be viewed as a 
transformation of the fundamental plane space. 

The use of L/2 as L1/2 would be invalid in the presence of 
significant attenuation due to dust, but the objects described 
here are have very low gas fractions and hence likely have 
negligible attenuation. Thus the interpretation of L1/2 = L{< 
ri/2) as the light emitted within ri/2 is a reasonable one for 
these objects, and the above set of observable transforms re- 
lations are close approximations to the actual MRL variables. 

Later, we will also consider a modified version of MRL 
space that we call dMRL space. In dMRL space, the mass 
variable is MP,^ = M1/2— -^'^baryon(< ri/2),ihs darkmatter 
mass within ri/2- For our purposes, the difference between 
il/R^ and Afi/2 will only be substantial for E and dE galax- 
ies, and is obtained by subtracting out the stellar mass within 
the half-light radius for these galaxies (which contain negli- 
gible gas fractions); expHcidy, MR^ ^ M1/2 - A/*/2. It is 
important to recognize that the presence of radial gradients in 
M^/L due to metallicity variation could render the use of our 
formula for My^ invalid by shifting the radius enclosing half 

the ste llar mass from ri/2- However, as shown in Smi th et al.l 
( 120091) . typical metallicity gradients for the E galaxies (for 
which Af* is most important) are 5 log{Z) / S log{r) « —0.1. 
Using this gradient with a typical anci ent (13.7 Gyr) solar 
metallicity stellar population from .Bruzual & Chariot, (,2003,), 
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we find M^jL shifts by 0.07 dex from i?e to 0.1_Re. Hence, 
this is a small effect for our galaxies and we disregard ifl- We 
return to dMRL space in the next section. 

In Figure|2]we plot the data set described in fJ2]transformed 
into MRL space and projected along the coordinate axes. For 
each of the projections, we have also plotted lines to reflect 
scalings of interest. 

The left panel represents £1/2 as a function of rxji- The 
galaxies show a trend of increasing luminosity with increas- 
ing ri/2' but occupy a relatively small fraction of the avail- 
able detection space. The GCs and UCDs, meanwhile, are 
much more scattered in this plot, and are consistently smaller 
than the dSphs at similar luminosities (e.g. higher surface 
brightness); as described below we interpret this (along with 
similar behavior in the other projections) as a clear sign they 
are separate populations. The dashed black line in the left 
panel is a line of constant surface brightness (ii/2 oc r\i^, 

fxv = 30 mag arcsec"^. Below this surface brightness limit, 
detection bias in this plane becomes signi ficant for MW dSphs 
jKoposov et all 12001 lWalshetal.1 12009). This likely biases 
the observed ?"i /2 ~ ^1/ 2 relation to small ri/2 at the faint 
end (iBullock et al.ll2010h . We discuss the effect of th is bi as 
on our parameterization of the 7-1/2 ^ ^1/2 relation in ^A.W 

In the middle panel we show a projection into the ri/2 — 
Mii2 space. We include lines of constant mass density 



(Ml/2 



1/2- 



black dash-dotted line) and constant surface 



mass density (M1/2 oc r^i2, black dotted line), with normal- 
izations arbitrarily set to bracket this data set. In almost all 
cases, spherical geometry is assumed, in which case a slope of 
2, is more properly characterized as a 3D density profile that 
varies as p oc r^^ (somewhat cuspier than constant density). 
A slope of 3, meanwhile, is the scaling expected if all galaxies 
had a single constant density within their half-li ght radii. This 
slope has b een noted previou s ly at some scales(|Gentile et al.l 
|I009; Nap ohtano et"al]l20Tol: IWalker et al.ll2010l) . The fact 
that the dSph galaxies lie above the constant density line 
(black dash-dotted) that is normalized to intersect the most 
massive cluster population suggests that they are slightly 
denser than galaxy clusters (but not that much denser) at their 
half-light radii. For a figure that explicitly compares the im- 
plied mean density of these objects, see Appendix IBJ 

Finally, in the right panel we show Afi/2 vs. £1/2, and 
a mass-follows-light line (Afi/2 oc £1/2) normalized at 
M/L = 3 in solar units to reflect the mass-to-light ratio of 
a uniform fairly old stellar population. Note that the deviation 
of a population from A/1/2 oc £1/2 is equivalent to the "tilt" 
that is often discussed in the context of the fundamental plane. 
It is clear from this figure that the CSphs and dSphs deviate 
from this scaling substantially owing to their high dark mat- 
ter fractions, while the other populations are more consistent, 
although the Es do show the well-known tilt, and the UCDs 
show a possible tilt (discussed below). 

Figure |3] shows the same information, now presented in a 
3-D r epresentation. The red plane outlined with a solid line 
is the Graves et al.l (|2009b) fundamental plane (transformed 
into MRL space). The blue plane outlined with a dashed 
line is a plane with mass proportional to luminosity with 
Ml/2 = 3£i/2 and is indicative of the plane we would expect 

^ In principle, a radial A/»/L .shift could be resolved by replacing Z/1/2 
by M^x/2 ^nd defining the appropriate ^,,1/2. However, the data quality is 
not sufficient to derive _A/, for our full sample, so we use L here. 



uniformly old, purely stellar systems to lie within. We note 
that, in fundamental plane space, this last scaling is sometimes 
called the "virial plane" (even though systems can be in virial 
equilibrium regardless of whether or not they lie within this 
plane). In MRL space it is manifestly apparent that this plane 
is defined by the assumption that mass-follows-light with a 
fixed M/L. 

Another feature revealed by examination of the populations 
in Figures |2]and|3]is a distinct separation between dSphs (ma- 
genta) along one sequence and UCDs/GCs (blue/green) along 
another (a similar situation is noted by Forbes et al. 20081 
in K-band). Specifically, the UCDs and GCs cluster more 
closely around the Afi/2 oc £1/2 plane (shown as dashed, 
transparent blue) while the dSphs (at similar luminosity) peel 
sharply up from it, reflecting a significant dark matter com- 
ponent and larger sizes. This difference is clearly visible in 
the two-dimensional projections of MRL space shown in Fig- 
ure 12] and manifests itself as a wishbone-shaped bifurcation 
of the spheroidal sequence in Figure [3] We also note here 
that the UCD sample seems to show a slight tilt from the 
A/1/2 oc £1/2 relation, most clearly apparent in the right 
panel of Figure |2] This could be a sign of a very small 
amount of dark matter, but could also be systematic varia- 
tion in the A/,/£ ratio due to stellar effects. These objects 
have uniquely large luminosity densities, and hence are the 
most likely places to show chang es in star formation condi- 
tions (" Dabringhausen et"an 120091) or simply be an extension 
of scalings that exist everywhere (such as variation for the Es 
is described in more detail in il4.3l l. Alternatively, they may be 
due to dy namical evo lution or more complex fo rmation sce- 
narios (e.g lGoerdtet a l. 2008; Taylor et al.l2010l) . Regardless, 
the significance of this tilt is not clear from this data set (al- 
though more tilted than the GCs), and the UCDs and GCs are 
quite distinct from the dSph sample. 

Given the observation that the M W dSphs ar e dark 
matter-dominated dSimon & Geha 2007; Strigariet al. 200813 
ISimon et al.l 120101) . and GCs have Af /£ consis tent with 
purely stellar systems (e.g. iPrvor & Mevlanlfr993h . we con- 
sider if there is a clean separation between these systems 
based on the MRL space parameters. We fit a plane that 
separates the dSphs from GCs by finding a plane that lies 
perpendicular to the best least squares fit between all of the 
dSphs and GCs and perpendicular to the best fit line through 
the dSph sequence; we then offset the plane until it evenly di- 
vides the two populations, giving the plane rendered in Figure 
ID This plane is a convenient empirical way to determine if an 
object is a faint dSph or a globular cluster In the MRL space 
for our data set, the best fit separation plane is given by 

0.34 log A/i/2-0.501og £1/2 + 0.79 log ri/2 = -1.35. (2) 

Specifically, objects that lie at lower A/1/2, lower ri/2, or 
higher £1/2 are GCs while others are galaxies. This same re- 
lation can easily be transformed into fundamental plane space, 
providing the separation plane 



0.68 log cr - 0.50 log /e + 0.13 log /?e 



-3.23. (3) 



such that objects with lower a, higher log I,,, or lower log /?e 
are GCs while others are galaxies. 

The fact that this single plane easily separates the GCs and 
dSphs in the MRL space implies that these are distinct classes 
of objects (see also the discussion in Appendix ICJ - the argu- 
ments there for UCDs also apply to GCs). It is possible that 
future studies of faint/low surface brightness GCs may change 




14 


- 


,!■■'#? 


^,12 


- 


R:f 


o 
^ 10 


# 


/■ 


^ „ 


/f> ,/ 




Tr » 


^■'•^ .' 


— 


o 






6 


'^ff^/ 


— 


4 


^ x'l* 1 1 •" 1 


1 - 



14 


- 


1 1 


' S 
f 


^ 12 


_ 




V 


^ 10 






y^ 


- 




/ 


^ 






^ 


^ 




f. 


p 










^ 




■^ 




6 




^ 


- 


4 


-•/ 


1 1 


1 



nan 


CSph 


OOo 


dE 


AAA 


dSph 


ft6<r 


GC 


• •• 


UCD 


• • . 


E 



-4-3-2-10 1 2 3 
log(ri/2/kpc) 



-4-3-2-10 1 2 3 
log(ri/2/kpc) 



4 6 8 10 

log(Li/2/LQ) 



12 



Fig. 2. — Two dimensional projections of the three-dimensional data shown in Figure[T]onto the coordinate axes of MRL space: luminosity vs. radius, mass 
vs. radius, and mass vs. luminosity (from left to right, respectively). The color/symbol code is indicated. In the left panel, the dashed line is a constant surface 
brightness of p^v = 30 mag/arcsec^ - any galaxy below this line would be undetectable with current methods. In the middle panel, the dash-dotted line shows 
the scaling of constant density (M oc r^) while the dotted line shows the scaling of constant surface density (M oc r-^), with normalization chosen to bracket 
the galactic objects. In the right panel, the sohd blue line reflects a constant mass-to-light ratio, specifically M1/2/MQ = 3LiI2/Lq. Slope deviations from 
the A/1/2 ^ ^1/2 li"2 ^2 equivalent to "tilts" in each galaxy population's fundamental plane. 



the location of this separation plane, or even fill in the gap, 
rendering the plane completely arbitrary. But for this data set, 
the classes are completely separated by the plane of FigureH) 
Further, we note that this plane implies that a galaxy/cluster 
projection using a single variable (e.g. iGilmore et al]l2007h 
is not sufficient to separate these populations, as is apparent 
from Figure m All 3 dimensions are necessary to account for 
the most extreme objects. 

Additionally, we include UCDs in Figure |4] and find that 
they also lie clearly separated by the plane, even though they 
are not included in the determination of the best-fit separa- 
tion plane. This is suggestive that they are in the same class 
as GCs, and not on the galaxy sequence. However, although 
given the tilt discussed above, we cannot discount the possi- 
bility that this is simply due to a relative rarity of the most 
massive UCDs to bridge the gap. 

Given that GCs and UCDs both lack clear evidence for 
dark matter and sit in a distinct region of MRL space we are 
inclined to treat them as stellar systems rather than "galax- 
ies", which we define operationally as stellar systems that are 
bound to a dominant dark matter halo (as discussed in S{T]i. 
Alternatively, a second scenario is possible where UCDs do 
contain significant dark matter If this is the case, then an 
interesting implication follows: there would need to be a di- 
chotomy in galaxy formation efficiency in dark matter halos 
of a fixed virial mass. Specifically, as shown in Appendix ICl 
most UCDs are consistent with no dark matter given the un- 
certainties in the expected stellar mass-to-light ratios. If we 
force a stellar mass-to-light ratio of 2 (such that their dark 
matter densities are comparable to their dynamical mass den- 
sities) then the implied dark matter densities are incredibly 
high - comparable to the central densities of the most mas- 
sive galaxy clusters (Mvir ~ IO^^Mq). dSphs of similar lu- 
minosities sit in i\/vir ~ IO^A/q halos. UCD dark matter 
mass fractions would need to be extremely fined-tuned (and 
different from object to object) in order to avoid a dichotomy 
in galaxy formation efficiency at a fixed dark matter halo mass 
- a dichotomy that is not seen for any other type of spheroidal 
system. This is an interesting possibility and may call for 
more investigation, as such a result would be difficult to ex- 
plain in LCDM. 

Nevertheless, we regard the above scenario to be unlikely, 
and adopt the simpler interpretation that UCDs are purely stel- 



lar systems that occasionally have unusually high M^,/L due 
to unique star formation conditions or dynamical evolution. 
From here on we omit the GCs and UCDs from consideration 
as systems that clearly contain dark matter halos of their own. 
In the alternative scenario where UCDs are to be regarded 
as galaxies, our approach could be viewed as restricting our- 
selves to the simpler dSph "branch" of the MRL relation. 

Once we remove UCDs and GCs, we are left with a galaxy 
sequence in Figures|2]and[3]that scatters about a 1-D relation 
through MRL space. In the next section we work towards 
characterizing this 1-D curve. 

Figure |5] provides yet another representation of the MRL 
data, now presented as the dynamical half-light mass-to-light 
ratio Ti/2 = ^iilil^xji (in Mq/ Lq) plotted as a function 
of each of the MRL variables individually. Along the top of 
each panel we show characteristic observational uncertainties 
for our galaxies of each type across the MRL sequence. We 
discuss these errors in the context of measuring scatter in the 
MRL relation in §6. 

Each panel in Figure |5] clearly reveals a minimum T1/2 — 
3 that spans a broad regime of spheroidal galaxies, from 



Ml/2 - lO^-iiMg (left); ri/2 c± 1 - 10 kpc (middle); 
andLi/2 ~ 1O''~^°L0 (right). As discussed by Wolfetal] 
(|20l3) in the context of a similar figure in their paper (Fig- 
ure 4), the dramatic increase in dynamical half-light mass- 
to-light ratios at both smaller and larger scales is likely in- 
dicative of a decrease in the efficiency of galaxy formation 
in the smallest and largest dark matter halos - as discussed 
above, the influence of radial variations in M^jL is ^ 0.1 dex, 
far less than that observed here. For the biggest, brightest, 
most massive galaxies, the increase in T1/2 implies a sharp 
threshold for galaxy formation in luminosity (not in mass) at 
Li/2 — 1O^^L0, as shown by the strong break in the left panel 
of Figure|5] The strong sensitivity to luminosity suggests that 
baryonic processes are responsible for this transition. Mean- 
while, the smallest, faintest, least massive galaxies seem to 
exhibit a sharp rise at a particular mass scale (not luminosity 
scale) near Ml/2 — lO^M© (right panel of Figure [S]). This 
indicates they are more tied to the size of their potential wells 
than star formation (although this does not preclude an inter- 
action between the two, e.g. Dekel & Wool 120031 ; IWoo et al] 
2008:) . We connect these scaling trends to dark matter halo 
virial masses and relate them broadly to galaxy formation in 
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Fig. 3. — Three dimensional representation of the MRL space of log(Mx m/A^©). log(ii/2/^Q). ^^d log(i?i /2/kpc). The transparent (red) plane with 
solid borders is the fundamental plane of lGraves et alj i2009b). The (blue) plane with dashed borders corresponds to A/j /2 = 3Li/2, i.e. the plane for mass- 
follows-light, sometimes refeired to as the "virial plane". The data point color and point-type scheme matches that of the Figure [T] (or see [|2). A rotating 
animation of this plot is available at,http: //www .physics .uci . edu/ -bul lock /f curve /mo vies .html 



Sections |5] and m 

4. FUNDAMENTAL CURVE 

It is evident in Figures |2] and [3] that CSphs, Es, dEs, 
and dSphs seem to curve through MRL space alon g a 1- 
D seq uence (see also Graham et al. 2006; Graham & Worleyl 
120081 for dE and Es). We refer to this sequence as the "fun- 
damental curve" and we plot analytic representations of this 
curve in the left panel of Figure |6| along with the associated 
data poi nts. We disc uss these analytic curve representations 
Sections l4TT]andl4!2l 

It is important to note that the existence of this 1-D 
curve does not imply that these objects are a single pa- 
rameter family, nor that the curve is a more suitable fit 
than a higher-dimensional construct. As the fundamental 
plane (Graves etal. 2009bi) for Es and fundamental mani- 
fold (" Zaritskv et al.i20 06b) show, galaxies do show systematic 
variation along multiple directions in fundamental plane or 



MRL space. We do not aim to compare the statistical signifi- 
cance of these relations to the fundamental curve, as the appli- 
cations of 1-D and 2-D relations are quite different. Instead, 
the best way to think of the fundamental curve is as the di- 
rection of largest variation of this set of dispersion-supported 
galaxy properties. Thus, it is useful as the first-order scaling 
relation, and hence the first priority is to understand galaxies' 
positions along the curve. The other significant scalings are 
then encoded in the "intrinsic scatt er" a bout the fundamental 
curve (discussed and quantified in i l4.3l and i l5.2l l. 

The right panel of Figure|6]shows the same data, but now in 
dMRL space. Recall that the only difference between dMRL 
space and MRL space is that the dynamical mass within the 
half-light radius, A/1/2, is replaced by the dark matter mass 
within the same radius: Afi/2 -^ ^'^Yl2- "^^^ half-light dark 
matter mass is determined by subtracting the stellar mass of 
each system via MYJ^ = Af 1/2 — M^^jl. For the E galaxies of 
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Fig. 4. — Three dimensional representation of the dSph/GC separation 
plane in MRL space. The green points are GCs, blue points are UCDs, and the 
yellow points are dSphs. The transparent cyan plane is given by Equation[2] 
chosen to separate the two populations, fit using only the dSphs and the GCs 
(se e text for details of the fit). A rotating animation of this plot is ava ilable at 
|http : //www .physics ■ uci . edu/ "bullock/f curve /movies . html| 



[Graves et al. ! ('2009b| ) we use stel lar masses derived from the 
estimates of Gallazzi et al.l (120051) (see^raves & Faber 2010, 
for more details). For the dE sample of Geha et al. (2003), 
explicitly computed stellar masses are unavailable so we as- 
sign them stellar masses f rom their o b served integrated colors 
using the prescription of iBell et alj (l2003h . For the CSphs 
and dSphs we assume Afi/2 = -^1/2'' because the dynamical 
mass-to-light ratios in these systems are very large. 

The motivation for exploring dMRL space and its funda- 
mental curve is that we would like to use the dark matter mass 
density within ri/2 as a estimator for the halo virial mass. 
With a virial mass estimate in hand, the fundamental curve 
relation can be used to provide an approximate, average rela- 
tionship between halo virial mass (A/vir) and galaxy luminos- 
ity (L). This necessitates comparison to a 1-D dMRL relation, 
as halo virial masses are a one-parameter family. We discuss 
this effort in E|5] 

4.1. MRL Curve Models 

We have chosen to quantify the fundamental curve by treat- 
ing ri/2 ^s the parametric variable. We fit two relations, one 
in the ri/2 — -^1/2 plane and another in the ri/2 ~ ^^1/2 plane. 
The derived pair of relations (RL and RM) define our funda- 
mental curve relation for the three MRL variables. We also fit 
the curve directly in three dimensions for some models, but 
the derived parameters were effectively identical, and hence 
we use the simpler two dimensional fits for clarity. We now 
describe our choice of functional forms for modeling these re- 
lations, followed by a set of five best-fit models for the funda- 
mental curve, distinguished by slight differences in the fitting 
procedure and the choice of M^J^ — M1/2 — M^jl as the 
mass variable in place of the raw Af^ /2. 



For the ri/2 ~ ^1/2 relation, we define f^ = log(ri/2/^L) 
and L = log(Li/2/-^o) and employ a fit following the 
empirically-motivated form 

L = rL^^ + [s-rL{a-b)] ^^^. (4) 



Equation |4] has the property of smoothly transitioning from 
an asymptotic slope a (such that ii/2 oc r^/2) for ^1/2 <^ r^ 

to 6 (i.e. Li/2 oc 7-^/2) ^'-'^ ''1/2 ^ ''"l, with the width of 
the transition zone at r^ defined by w. The parameter Lq 
is then the characteristic luminosity at r = r^,, and the final 
parameter s determines the size of a luminosity offset that 
occurs in the transition region (e.g. the break in luminosity at 
log(ri/2) ~ in the upper-middle panel of Figure |7]l. This 
fitting function simply yet generic ally captures the behavior 
of a data set that has distinct asymptotic power laws and a 
smooth transition region between them. 

For the ri/2 — M1/2 relation we utilize a fitting function 
with a form identical to Equation|4l 

^ a. + j3 , . , ^,, arctan(fM/i^) 

Z TT 

where tm = log(ri/2/?'Af ), so that r^ defines the transition 
radius and A/ = log(Afi/2/A^o) defines a characteristic mass 
scale Mq at r — tm- 

Using this method, the A/1/2 vs. L1/2 relations are gen- 
erated by eliminating our chosen parametric variable ri/2 in 
Equations |4] and |5] For comparison, we also directly fit the 
ML relation using the form of EquationH] and find very sim- 
ilar relations to those shown below. Hence, the results pre- 
sented here are likely not very sensitive to the choice of ri/2 
as the parametric variable. 

Motivated by the fact that we are interested in understand- 
ing each type of galaxy universally (CSph, E, dE, and dSph) 
we weight the data points such that each of the four groups 
has equal weight (i.e. the weight for each point is 1/Ntype 
where Nttjpe are the number of objects of that type). Further- 
more, for the E data set of Graves et al. (2009b), we weight 
each point by the relative fraction of galaxies in that partic- 
ular bin so as to properly represent the full SDSS population 
rather than the choice of bin locations. With these weights for 
the data set, a non-linear least-squared fit for the parameters 
in Equations|4]and|5](using a LevenbergMarquardt algorithm) 
fully determines the one-dimensional relations. 

With the models for the fundamental curve and this fit- 
ting procedure, we call our empirically fit fundamental curve 
model "MRL- 1 ," with best-fit parameters given in the first col- 
umn of Table |2] The relation is shown in projection on the 
MRL axes as a blue-dashed line on the top panels of Figure 
|7] along with data points for the individual galaxies and their 
associate d ob servational error bars (error bars are discussed in 
detail in il4.3l l. We show the same curves and data points as 3- 
D representations in|6] with the MRL-1 shown as the dashed 
green line in the left panel. 

The black dotted line in the upper left panel of Figure 
I2] shows the surface brightness detection limit for dSphs, 
/i^^= 30 mag arcsec^^ (Koposov et al. 2008; Walsh et aQ 
l2009h . Given that the detection limit indicates that the least 
luminous dSphs galaxies are at the edge of detectability, it 
is plausible that the shallow slope in RL at fai nt £1/2 is due 
to a selection effect. The "stealth galaxies" of iBullock et alj 
(,201(X) . if present, could substantially alter the slope at the 
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Fig. 5. — Mass-to-light ratios witliin tlie lialf-light radius, Tj /2 = Afj /2/L1 /2 in units of Mq/Lq, shown as a function of each of the MRL variables. 
Error bars shown along the top of each panel are representative of the observational uncertainties in each parameter and each galaxy type, indicated by matching 
color-code and location in the M, R, or L, axis. For the dSphs, we include separate error bars for the classical and SDSS dwarfs. The "?" for the CSph M1/2 
indicates the additional (unquantifiable) uncertainty unique to the CSph, due to the use of the cluster galaxies to detennine the velocity dispersion instead of the 
actual dispersion of the ICS. 



faint end. Thus, we also include an "MRL-2" model in which 
the s parameter is forced to be 0, causing the faint end slope to 
trace the full dSph population instead of being strongly driven 
by the faintest of them. This model is shown as the solid black 
line on the left panel of Figure|6]and the upper panels of]?] and 
the best-fit parameters are given in the second column of Ta- 
ble|2] Given the fact that r nost of the faint dSphs skirt the edge 
of this detection limit (e.g. lWalsh et al.l2b09,) . we consider the 
MRL-2 model to be the more robust choice for characterizing 
the MRL fundamental curve. 

The fit parameters listed in Table |2] for the MRL-2 model 
reveal that the smallest galaxies with L < 2 Lq ~ 4 x 10^ Lq 
follow a mass-luminosity relationship that varies weakly with 
luminosity 



Ml/2 oc L 



a/a ^ 7-0.30 
1/2 ^ -^1/2 ■ 



(6) 



while the largest galaxies (i > 4 
mass-luminosity relationship with 



10^ i©) obey a steep 



Ml/; 



OC i^y2 °^ -^1/2 • 



(7) 



Both regimes are clearly very far from mass-foUows-Ught 
scalingfl(i.e., M1/2 oc L1/2). 

For the smallest galaxies, large changes in luminosity corre- 
spond to fairly minor changes in half-light mass. Conversely, 
for the largest galaxies, a factor of ~ 2 change in luminos- 
ity corresponds to more than an order of magnitude change 
in half-light mass. This is the same effect noted in Sj3](with 
regard to Figure|5]l, and without any appeal to theory suggests 
that two qualitatively different processes are acting to sup- 
press baryon conversion into stars along the transition from 
small galaxies to large. The smallest galaxies seem to be lim- 
ited by the dark matter mass itself (e.g., by the potential well 
depth), while the largest galaxies seem to be baryon limited 
(e.g., by the supply of cool gas for star formation). 

Also of interest is the sharp transition in the RM relation 
at log(ri/2) — 0.5 and \og{Mi/2) — 9, where the half-light 
mass suddenly jumps with increasing radius. This transition 
scale corresponds closely to the point where the dynamical 
mass-to-light ratios of galaxies reach their minimum (Figure 
ID and thus where baryons contribute substantially to the mass 

* For r^ and rjv/ values somewhat different from the best-fit for this data 
set, the values of these slopes can be quite different, but mass-follows-light 
never holds for any reasonable fits. 



compared to dark matter. It is possible that this feature is en- 
hanced or even caused b y the effects of baryonic contraction 
jBlumenthal et al.lll986l) as discussed in the context of dark 
matter masses below. 

4.2. dMRL Curve Models 

Recall that the dMRL relation is distinguished from the 
MRL relation by the use of M^J^ ^ A/1/2 - M^/2 as the 
mass variable in place of the raw dynamical A'/i/2. The fit to 
the data in this space using Equations |4] and |5] is our "dMRL- 
1" model. Trying a variety of starting values for the parame- 
ters revealed that rj^i is not well-constrained by the data and 
often would end up outside the data set regardless of the start- 
ing value. Hence, we used the RL relation to set the scale, 
through the constraint tm — r^. Using this constraint, the 
final parameters are given in the third column of Table |2] and 
plotted in the right panel of Figure |6] and the lower panels of 
Figure|7]as the red dashed line. 

As Table 12] shows, the dMRL-1 model best- fit parameters 
have a K, 13, and a preferring 0. Equation |5] for dMRL-1 
reduces to a power law for a — (3 and cr = 0, so the ri/2 — 
Af P.2' relation turns out to be very close to a single power law 

(linear in M and tm)- Hence, the RM relation can be modeled 
as a simple power law 



AfDM 
^*''l/2 



A/o(ri/2/rM)", 



(8) 



where r7\/ is determined from the dMRL- 1 fit to simplify com- 
parisons. The value of the slope a — 2.33 is also given in 
Table|2] The lower-middle panel of Figure|7]compares this fit 
(red dotted line) to dMRL- 1 , showing an insignificant differ- 
ence. 

Thus in the second dMRL model (dMRL-2) we adopt Equa- 
tion [U as the model for the RM relation, and the RL model 
of MRL-2, selected due to the likely presence of the stealth 
galaxy selection effect. We tabulate the best-fit parameters 
for this model in the second-to-last column of Table |2] and 
plot it as the black solid line in the lower panels of Figure [T] 
and the right panel of Figure |6] 

In the RM re lation of the dMRL space, we include for 
comparison the IWaUcer et al.l (12010 ) relation derived using 
Milky Way dSphs for the faint end and spiral galaxy rota- 
tion curves for the galaxy regime (black dash-dotted line on 
the lower-middle panel of Figure |2]l. We note here that while 
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Fig. 6. — Three dimensional representations of the best-fit fundamental curve relations in MRL space (left panel) and dMRL space (right panel, with Mj /2 — > 

J\/p,^). Our fiducial models MRL-2 (left panel) and dMRL-2 (right panel) are shown by black solid Unes. Models MRL-1 (left panel) and dMRL-3 (right panel) 

are shown by a dotted green line. The data point color and point-type scheme matches that of t he previous figures (see ij2j. Rotating animations of these plots are 
available at|http : //www .physics . uci . edu /~bul lock /f curve /movies . html] 
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Fig. 7. — Projections of galaxy data and fundamental curve in MRL (upper panels) and dMRL (lower & upper-left panels) spaces. The best-fit fundamental 
curve models for MRL-1 (upper, green dashed line), MRL-2 (upper, solid black line, fiducial) dMRL-3(lower, green dashed fine) and dMRL-2 (lower, sohd black 
line, fiducial) are shown in each projection. The straight dotted line in the left-hand panel represents the detection completeness limit for dSph galaxies found in 
the SP SS of {lY = 30 mag arcsec~^. The black dash-dotted line in th e lower middle panel is the relation for for dSphs and spiral galaxies from Walker et al] 
1201(1) . Enor bars are observational error bars for each galaxy (see i|4.3l for more details). The data point color and point-type scheme matches that of the previous 
figures (see ij2). 
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the IWalker et alj (120101) non-dSph sample is a very different 
set of galaxies that may obey different scaling relations from 
our sampl^E it is fairly close to our relation in the galaxy 
regime. However, the relation steepens with the inclusion of 
Es and CSphs, so our derived slope is somewhat higher than 



TABLE 2 

Fundamental Curve Model Parameter Values for Equations 

|4]and[5] 



a Af °^ oc r^ relation. 
Motivated partly by this M 



DM 



r-^ result on the faint end. 



as well as the greater uncertainty in M^,^ for the dEs and Es 
(see g3]and|52]i, we consider a third dMRL model (dMRL- 
3). In this model we use Equation |5] for the RM relation, but 
we force the faint-end slope to 2 and set the normalization to 
pass through the dSphs. We then force the rj\/ scale to match 
r^ (from MRL-2), set co = 0.01 to ensure a small transition 
rgion, and fir the remaining parameters. We also continue to 
use the RL model of MRL-2 for dMRL-3. In the last column 
of Table 12] we show the best-fit parameters of this model, and 
in the lower panels of Figure |7] and the right panel of Figure 
|6] we plot it as the green dashed line. 

Before continuing, we note a discrepancy for the E galax- 
ies in the dMRL models, most apparent in the lower-middle 
panel of Figure |7]- the Es tend to have higher A/-R^ than the 
best-fit relations. Recall, however, that the primary motivation 
for exploring the A/P.^ as a parameter is that it will allow us 
to map galaxy properties to an underlying dark matter halo 
mass. This mapping is hindered somewhat by the contrac- 
tion of baryons. An anomalously high dark matter mass for 
the galaxies with the highest baryonic-to-dark matter ratio is 
precisely what is expected if dark matter halos contract due 
to cen tral condensation of baryonic matter (iBlumenthal et al.l 
Il986h . Thus, we might expect an offset in the scaling relations 
of galaxies at the scale where baryonic condensation has been 
the most significant. In i|5]we estimate the degree to which 
baryonic contraction has increased the M^^ masses in our E 
galaxy sample and show that this increase approximately ac- 
counts for the discrepancy. Further, as discussed more in il4.31 
a power law for the RM relation is in general more robust to 
the problem of a non-monotonic mapping of baryonic galax- 
ies to dark matter halos. Thus, use of a power law for the RM 
model is a reasonable choice for the exercise of halo profile 
matching (described in fJ5]l, while still being a decent fit to this 
data set. In the RL space, as described above for MRL-2, it is 
more appropriate to use the s ~ model so as to prevent the 
stealth galaxies selection effect from strongly biasing the faint 
end slope. Thus, we adopt dMRL-2 as our fiducial model in 
the latter sections of this paper. 

4.3. Scatter and Uncertainty in the Fundamental Curve 

It is interesting to ask about the degree of intrinsic scatter 
within the fundamental curve that was defined in the previous 
section, but in order to do that we need to estimate the obser- 
vational uncertainties on the MRL variables. Representative 
error bars for M1/2, ri/2, and L1/2 are shown in Figure|5]for 
several galaxy types. Observational errors for M^J^ are pre- 
sented in Figure |7] Individual error bars for each data point 
are shown in Figure [T] Note that for the faint dSphs and the 
CSph, the measured mass-to-light ratios are much larger than 
any reasonable stellar population (e.g. M1/2/L1/2 » !)■ 
Hence, they are dark matter-dominated {M1/2 ~ AfR^), and 
hence the M1/2 and M^Jf errors are similar to each other. 

' See lMcGaugh & Woljj2010h for a discussion of liow dSpli scaling rela- 
tions connect to spirals. 
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For the dE and E galaxies, however, the mass-to-light ratios 
are closer to that expected of stellar populations and hence a 
significant amount of mass within ri/2 is in stars rather than 
dark matter, so M-R^ errors are larger for these objects due to 

the errors on Af * . 

For the E galaxies, the uncertainty in M^ due to stellar 
populations is a major uncertainty. While the observational 
errors play a role in general, for the large stacked E data 
sets here, the errors are certainly dominated by systematics, 
of wh ich there are three major components dGraves & Fabeij 
I2OIOI) . First, there is variation due to the method used to de- 
rive M^jL (e.g. int egrated colors or par ticular spectral fea- 
tures). As shown in [Graves & FabeJ (|2010), this contributes 
a Icr scatter of ~ 0.08 dex. Second, the assumed star for- 
mation history affects the inferred stellar mass, a t the level of 
0.15 dex for this sample (iGraves & Fabeiil2010l) . Third, the 
choice of IMF has a major effect on the inferred M*. For 
the example (conservative) comparison of Chabrier as com- 
pared to Kroupa (Longhetti & Saracco 2009), the inferred M^ 
varies by 0.26 dex. More detailed studies of individual objects 
can p otentially reduce the systematics (e.g. Cappellari et alj 
I2OO6I) . but the analysis above is appropriate for the large data 
set in use here. Thus we show error bars by adding the above 3 
components in quadrature, providing a factor of 2 uncertainty 
in the Af* used for mapping Afi/2 to AfP,^. This error on the 
Es is shown in Figure|2]as the error bar on A/R^, and we also 
adopt it in the next sections as the error for A/P.^. 

The error bars shown in Figure|2]account for the uncertainty 
in measuring the dark matter mass as it is today, but do not in- 
clude the systematic uncertainty that remains in our ability to 
map an observed dark matter density to the virial properties of 
that d ark matter halo. Baryonic contraction ( Blumenthal et al] 
119861) in particular can make the mapping between density 
and global halo mass quite difficult. We expect this uncer- 
tainty to be particularly important for E galaxies because they 
have the highest baryon fractions. We discuss this effect in 
more detail in iJS] 

For the dE sample. A/* is inferred from SDSS colors as de- 
scribed in 121 Errors can be estimated from this procedure by 
comparing the inferred Af* for each band. Using this proce- 
dure, the scatter in the inferred A/* is about 30%, comparable 
to the observational errors for Afi/2- This estimate has its own 
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set of systematic errors like those described above - we do not 
quantify this here due to the smaller sample size (and hence 
larger random errors) and more simplistic method compared 
to the Es. Regardless, the error bars are large enough to be 
consistent with the fundamental curve. 

For the CSph population, the uncertainty in Afi/2 is diffi- 
cult to characterize, as it is primarily due to the use of the 
galaxies to trace the velocity dispersion instead of the ICS. 
The effect this will have is not as well understood, as repre- 
sented by the "?" in the error bar of Figure|5] The simulations 
of Dolag et al. (2010) find a disagreement in cr of ^ 20% be- 
tween galaxies and the ICS component (i.e. approximately 
50% in mass), although this is not necessarily representative 
of the clusters in our sample. In order to broadly represent 
this uncertainty, we have assumed a factor of 2 uncertainty on 
G in deriving the error bars in the next section. 

Adopting these observational error bars, it is clear from the 
top panel of Figure |7] that the actual scatter about the fun- 
damental curve is larger than the observational errors. We 
estimate the scatter by computing the residuals of M^ji and 
^1/2 from the fundamental curve, and measure the standard 
deviation with weights as described in i|4] The resulting as- 
observed scatter in A/1/2 at fixed ri/2 about the MRL-2 re- 
lation is (5 log Ml/2 = 0.41. Subtracting the observational 
error in -/\/i/2 (including the contribution due to the error 
in ^'1/2) in quadrature from this value, we obtain an esti- 
mated intrinsic scatter of AlogAfi/2 = 0.38. Using the 
dMRL-2 relation, the observed scatter in A/P,^ at fixed ri/2 is 
(SlogAff/^ = 0.60 and the intrinsic is AlogA^^^^ = 0.20, 
due to the larger uncertainties in MYj^ ■ We emphasize that 
this estimate of intrinsic scatter is only approximate, given 
our small samples size and our rough characterization of ob- 
servational errors over the entire (disjoint) population of our 
objects. 

Applying the same method to the RL relation (identical for 
MRL-2 and dMRL-2), we get an observed scatter in _Li/2 at 
fixed ri/2 of (51ogii/2 — 0.73, and estimated intrinsic scat- 
ter AlogLi/2 = 0.71. This is relatively high, but is driven 
almost entirely by a few dSph outliers (the low dSph points 
in the upper-left panel of Figure |7]i that render the distribu- 
tion non-Gaussian. The dSphs generally have relatively high 
error bars, but this is not accounted for in the averaging pro- 
cess above. Thus, removing the discrepant dSphs gives an 
observed scatter of L1/2 at fixed ri/2 of (51ogLi/2 = 0.42, 
and A log L1/2 = 0.37. 

These values for the scatter are purely empirical measure- 
ments of the deviation of individual galaxies from the funda- 
mental curve. As discussed in fJ4] the intrinsic portion of this 
scatter encodes all of the additional scalings in galaxy for- 
mation that are sub-dominant to the curve itself. In the next 
section, we describe theoretically expected scatter based on 
the profile matching scheme. 

5. DARK MATTER PROFILE MATCHING 

We now describe a technique to use the fundamental curve 
described in the last section to derive global relations connect- 
ing dark matter halos to the luminous properties of the galaxy. 
The main relationship we would like to derive is the median 
relation between Afvii and L. We refer to this method as "pro- 
file matching," as it matches the mass profile of galaxies to 
dark matter halos to do t his. While the analysis presen ted 
here relies on NFW halos dNavarro. Frenk. & Whitell 19971) in 
ACDM, the general approach is applicable to any halo form 



or variant cosmology. 

ACDM simulations predict that at a fixed physical radius r, 
a more massive dark matter halo will be denser, on average, 
than a less massive dark matter halo (e.g. Navarro et al. 1997). 
Moreover, the typical mass profile for a given virial mass halo 
is determined by the virial mass in a one-to-one way, such that 
knowledge of AfP.^ and ri/2 for a galaxy can be mapped to 
the unique dark matter halo virial mass that gives Afhaio('' = 
ri/2; Afvir) = -^-^i/2^- Of course, this mapping is not without 
scatter, and we address this issue in ij5.2| This mapping is also 
made more difficult by the fact that some of the galaxies we 
consider reside within subhalos. We also address this point in 



We assume that each galaxy resides at the center of a dark 
matter halo and that galaxies have AfP,^, ri/2, and L^ii val- 
ues specified by the dMRL fundamental curve. We also as- 
sume that the dark matter densities within ri/2 can be mapped 
to a virial mass using density scaling relations derived for dark 
matter halos from dissipationless simulations. This is a rea- 
sonable assumption for most of our galaxies because most of 
them are dark matter dominated. This is not a good assump- 
tion for E galaxies, which have fairly high baryon mass frac- 
tions and have likely had their dark matter masses enhanced 
within ri /2 by baryonic conti'action dBlumenthal et alj[l986l : 
iGnedin et alJl2004lNapoUtano et al.ll2010l) . But as discussed 
in the previous section, the dMRL curves tend to lie below the 
dark matter masses in E galaxies in dMRL space. Indeed, we 
will show that a first-order correction for the effects of bary- 
onic contraction yields "uncontracted" masses for E galaxies 
that sit along our dMRL fits. 

We consider an ensemble of dark matter halos with a range 
of virial masses 10^ < Afvir/A/© < lO^*^-^. Each halo 
is assumed to follow an NFW mass profile Afhaio('') = 
Af (?■; A/vir) with a concentration parameter (c = rvii/rs) 
set by the median concentration-mass relations provided by 
iKlypin et al.l (120101) from the Bolshoi simulations. This sim- 
ulation was run with cosmological parameters {n — 0.95, 
0-8 = 0.82, h = 0.7, and 0,„ = 0.27) that are ver y 
similai- to those favored by WMAP7 dKomatsu et al.ll2010l). 
We de fine virial mass and virial radius as in Klvpin et alj 
(120101) . using the virial overdensity as calculated by the 
spherical collapse approximation. Note that we have ex- 
trapolated their fitted concentration-mass-redshift relation to 
masses beyond those directly probed by the Bolshoi simu- 
lation (A/vir = 10^"-^^^*-^Afo). However, these extrapola- 
tions are consistent with the scaling behaviors expected from 
previous simulations t hat have probed higher and lower mass 
regime s directly (e .g. iNeto et aljr2007l: ISpringeletaP 120081: 
MacciCEaDHOol). 

The implied dark matter mass profiles for many different 
virial masses are illustrated as colored lines in Figure [8] For 
reference, the half-mass radii for the dark matter halos, i?i/2°, 
are plotted as large colored circles at their associated half- 
mass values, Af j^f2° — A/vir/2. The slope of this relation is 
almost exactly Af j^?2° °c (i?5'/2°)^i and therefore significantly 
steeper than the M^J^ oc r\-h slope favored by our fiducial fit 
to the fundamental curve of stellar systems. The virial mass 
associated with each mass profile plotted is indicated to the 
right of the associated colored circle. 

Overlaid on Figure [8] as a thick, black solid line is the 
M^J^ vs. ri/2 relation for our preferred fundamental curve 
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Fig. 8. — Profile matching schematic. The solid colored lines are the enclosed mass as a function of distance from the center for NFW profiles of a given mass 
assuming a c-AI relation as described in the text. Their virial radii are enumerated at R^j}? (i.e. M(K!?^'°) = Mvir/2) as large colored points. The black solid 
line is our fiducial fundamental curve relation (dMRL-3) projected into this space, while the green dashed line is dMRL-3 (see [|4). The points are the data points 
used to fit the fundamental curve (see ^for meaning of colors and symbols). Observational error bars are shown for MP,^ and r^/g, derived as described in 
JI4.3I We note that the observational error bars on each point are significant and this contributes to the apparent scatter in the data. For the E data points, we 

rDM 
-'1/2 



correct the observed MPJ^ for baryonic contraction using the lBlumenthal et arHl986l) adiabatic contraction estimate described in the text, which results in an 



offset of 

the corrected values. 



5 dex in log(Afp,2 )■ In ths inset, we show the E galaxies before this correction as red open circles, and the solid red points with black outhnes are 
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Fig. 9. — Global properties of dispersion-supported galaxies embedded in 

dark matter halos. All panels have virial mass on the horizontal axes. The 

vertical axes are half-light mass (upper-left), half-luminosity (upper-right), 

3-D half-light radius (lower-left), and virial mass-to-light ratio. The blue 

thick lines with points are the inferred relations for the dM RL- 2 model, and 

the dotted red hnes are the associated fits using Equation lAll and Table \3\ 

The thick green dashed line corresponds to the dMRL-3 model. The black 

dash-dotted line in the upper-left panel shows Afj /2 = /baryon^vir- The 

1/3 
black dash-dotted line in the lower-left panel shows rj /2 oc M^/^. , which is 

the scaling one might naively expect in the case that stellar radii scale self- 

similarly with virial radii ri/2 oc iJvir- 

fit (Model dMRL-2 in Table IS . The thick green, dashed line 
is the alternative dMRL-3 relation. Each point along these 
curves is mapped to a single luminosity via its respective 
dMRL relation. Each point on the line can also be mapped 
in a one-to-one way to a median dark matter halo virial mass, 
set by the particular Afhaio(^) — M{r; Afvir) halo line it in- 
tersects. This allows us to back out an implied median re- 
lationship between galaxy luminosity and halo virial masses 
across the range of galaxies considered. Figure |9] shows the 
implied A/vir — L mapping for each of these curves (dMRL- 
2, solid blue with points and dMRL-3, green dashed) in the 
upper right panel. Associated relationships between A/vir and 
the other fundamental curve parameters are shown in the other 
panels of Figure |9] Full analytic descriptions of these rela- 
tion are provided in Appendix lAl (see Table O. For dMRL- 
2, the Ml/2 vs. A/vir and ri/2 vs, A/vir relations are fairly 
well characterized by power-laws with A/1/2 c« M^^^^ and 
ri/2 oc M^^^^. The L-to-A/vir relation, meanwhile, can be 
approximated on the faint end as i oc M^^^'^ and on the bright 
end as L oc M^^f^. As expected from our A/1/2 ^ ^1/2 deal- 
ings, one interpretation is that mass is the limiting factor in 
galaxy formation for faint galaxies while baryonic feedback 
of some kind limits galaxy formation for bright galaxies. 

Returning to Figure |8] we have also plotted the galaxy data 
points used to fit the fundamental curve as colored symbols, 
with error bars reproduced from the lower middle panel of 
Figure|2] The symbol types are the same as those described in 
ij2]and Figures |5]|2]except for the red (E) points, as described 
below. Clearly, these points exhibit a large scatter at fixed 



radius. As we discuss (and illustrate) in the next section, one 
of the reasons for the apparent scatter and offsets is that the 
measurement errors on each data point are quite large. This 
is particularly important for the red symbols (Es), for which 
small errors in stellar mass estimation can propagate to very 
large errors in the dark matter masses plotted, potentially in a 
systematic way. We discuss inherent vs observational scatter 
in detail in i j5.2l 

Another effect that adds uncertainty to the mapping be- 
tween halo mass and galaxy luminosity is baryonic contrac- 
tion (Blumenthal et al. 1986; Gnedinetal. 2004), which in- 
creases the dark matter density within a given radius from 
what it otherwise would have been absent the infall of 
baryons. The E points (red circles) in Figure|8]have been mod- 
ified in their A/P,^ masses from those shown in Figures|6]and 
I2] in order to approximately account for this effect. Specifi- 
cally, the DM masses for the E galaxies in this plot are esti- 
mates of the "intrinsic" dark matter masses within ri/2 prior 
to the infall of baryons. We mak e this estimate using the con- 
tra code of Gnedin et al. (2004) applied to the E galaxy bin 
with the largest number of galaxies. 

In order estimate the degree of the mass enhancement from 
baryonic contraction, we assume that the initial virial mass 
followed is that implied by our fiducial curve in Figure |9] 
(dMRL-2) for the ri/2 of the chosen E bin. We use the 
concentration-mass relation discussed above to det ermine the 
Cvir for an NEW profile. For simplicity we assume a lHernquisl 
(1990) model for the stellar distribution with A/* and ri/2 
set by the E bin. We determine the ratio of the mass within 
ri/2 before and after the contraction, and correct our profile 
matching A/R^ by this ratio. The points shown in Figure 

[8] assume the iBlumenthal et al.l ( 119861) adiabatic contraction 
formula, bu t we find t hat w ith both the lGnedin et al.l (120041) 
and Blumen thal et al.l (Il986h methods, the correction is large 
enough to move the E galaxies onto the dMRL-2 relation. 
For simplicity, the error bars on the E points here are simply 
scaled versions of the direct uncertainty in MR^ as presented 
in Figure |7] and do not include the additional uncertainty in 
the baryonic contraction correction, which is certainly large 
but hard to quantify. The errors shown here are conservatively 
small for this reason. 

The uncertainty in profile matching in the E/dE regime is 
nicely illustrated by the differences between the solid curve 
(from dMRL-2) and green dashed curves (from dMRL-3) in 
Figure |9] The dMRL-3 relation yields bumps (e.g. a plateau 
in L around A/vir ~ 10^^ A/© ) due to the enhanced MR^ at 

log(ri/2) ^ associated with this relation. This break in the 
MR relation maps onto an increased A/vir, creating this unex- 
pected feature, which is likely an artifact of baryonic contrac- 
tion, possibly with a component due to uncertainties in A/,. 

Regardless of the nature of this bump, however, this dMRL- 
3 scaling does a slightly better job in matching the properties 
of the faintest galaxies, as it was designed to have an MR rela- 
tion that is overweighted in dSph regime (compare the dashed 
and solid lines in Figure [8]l. Interestingly, the green dashed 
curves in Figure |9] reveal features in the scaling relations of 
the smallest galaxies at A/vir ~ lO^Af0 in the form of a 
wall in A/vir- Strictly speaking, this is a breakdown in mono- 
tonicity of the L — A/vir relation (discussed further in ^S.li . 
but for dMRL-3 this is because A/vir is very nearly constant 
with L. This m ight be indicative of a common mass scale for 
small galaxies (IStrigari et al.ll2008at iPenarrubia et alJl2008at 
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lOkamoto & Frenkll2009t:rWolf et al.ll20Tol) under which lumi- 
nous galaxies do not inhabit dark matter halos. Abundance 
matching does not constrain the existence of such a scale, as 
the galaxies in those halos are too faint to be observed in sta- 
tistically significant quantities outside the Local Group. As 
we discuss below, profile matching is just approaching the 
point where we can begin to test this possibility as part of 
a global relation. 

5.1. Comparison to Abundance Matching 

Figure [TO] compares our fiducial profile matched results 
(blue lines, dMRL-2) to those of the independent technique of 
abundance matching (red lines). The implied ratios (A/vir — 
L) vs. Mvii are shown in the left panel and the equiva- 
lent relations for (Afvii — L) vs. L are shown in the right 
panel. The blue profiled-matching lines are shown as dashed 
in the regime where the average dynamical mass-to-light ra- 
tio within ri/2 is indicative of a significant stellar component, 
with M1/2/L1/2 < 9. The line is solid in the regime where 
our stellar mass subtraction is less important for the dark mat- 
ter mass determination within ri/2- The line types emphasize 
the point that our profile matching technique is most trustwor- 
thy in the luminosity/mass extremes. We return to this point 
again in i i5.2l 

The red curves, specifically, illustrate the A/vir — L relation 
that is set by forcing the cumulative abundance of dark matter 
halos more massive than A/vii to match the observed cumula- 
tive abundance of all galaxies brighter than L (described, for 
example, in Kravtsov et al. 2004; Conrov & Wechsler 2008; 
iBusha et all I2009t |Moster et al. 2009). ' We use the SDSS 
luminosity func tion of Blanton et al. (2005) and the halo 
mass function of [Tinker et al.l (120081 for WMAP7 cosmolog- 
ical paramete r s). To convert from the SDSS bands used in 
iBlanton et al.l (120051) to the F-band used in this work, we 
use the transformation V — g — 0.59 * [g — r) — 0.01 from 
lister et al. (2005), implicitly assuming all galaxies have av- 
erage colors. The line becomes dashed where we have extrap- 
olated beyond the luminosity function completeness limit and 
becomes dotted at large luminosities where statistical uncer- 
tainties affect our ability to quantify the luminosity function. 

It is encouraging in Figure [T0| that our derived profile 
matching relation for dMRL-2 (blue, with circles) reveals 
a similar U-shape as does the abundance matching relation 
(red). In particular, our profile matched curve reveals a 
minimum of M^^^/L ~ 80 at Afvir — 2 x lO^^Af© and 
L 2± 2 X 10^°Lo, reflecting scales where galaxy formation 
efficiency is maximized. Similarly, the abundance-matched 
curve minimizes at AI^^^/L ~ 80 at A/vir — 3 x lO^^A/0 
and i ~ 4 X IQ'^Lq. This factor of ~ 6 agreement is rea- 
sonably encouraging, considering that the minimization of 
the abundance-matched curve occurs well within the regime 
where abundance matching is most affected by baryonic un- 
certainties. Compare the minima to the mass-to-light ra- 
tio that would result in the limiting case where 100% of 
each halo's baryons is converted to stars: (A/vir/i)inin — 
T*//baryon — 12 with T, w 2 Set by the average stellar 
mass of the E sample in this work (T^ = 1.89). The range 
1 < T, < 3 is shown in Figure[TO]as the gray shaded region 
clearly below any of the matching curves. The implication 
is that even for galaxies that are maximally efficient in con- 
verting their baryons into stars, some ~ 70% of their baryons 
remain unconverted. Of course, the inefficiency of baryon 
conversion into stars is a well-known result of CDM-based 



comparisons to galaxy luminosity functions. Nevertheless, it 
is encouraging that our profile matching analysis seems to im- 
ply the same level of inefficiency (on average) without appeal- 
ing to abundance information in any way. 

While the broad-brush agreement between abundance 
matching and profile matching is encouraging, clearly dis- 
tinct differences are present for dMRL-2. There could be 
several explanations for this. The most straightforward is 
that our profile matching M^-^jL relations are applicable to 
dispersion-supported galaxies, while abundance matching ap- 
plies to galaxies of all types. This is particularly important in 
the mass range A/vir — 10^°^^'^ A/© where the population of 
disky late-type galaxies become much more important relative 
to spheroidal early-types as mass decreases. The star form- 
ing galaxies will have higher luminosities (lower A/vir //^) 
than their pressure-supported/passive counterparts at the same 
A/vir, and it is only this latter category that is reflected in our 
profile matching data set. Hence, if the star formation effi- 
ciency peaks at a different mass for early-type galaxies than 
late-types, the two methods will give different results for the 
galaxies in this mass range. 

Additionally, at the bright end, abundance matching typi- 
cally matches the largest dark matter halos to bright E galax- 
ies. Thus they do not include the more diffuse, harder to mea- 
sure intra-cluster stars. We have included the full CSph light, 
and therefore the profile matched relation has a larger L at 
fixed A/vir (or lower My-^jL). 

With this in mind, it is important to note that at the cluster 
scale, direct object-by-object comparisons of the measured ef- 
ficiency (G onzalez et al. 2007) is complimentary to the scal- 
ing relation approach for comparison to galaxy formation 
models. Further, it is possible to directly co mpare lensing- 
based mass estimates to the stellar mass (e.g. IZaritskv et al.l 
2008). With a large enough sample, this could potentially de- 
termine whether there is a discrepancy in either abundance 
matching or profile matching, although the abundance match- 
ing estimates are rather uncertain at these mass ranges due 
to the impact of small numbers of large clusters (discussed 
above). However, because clusters are, by nature, systems 
where the subhalos/lower-luminosity galaxies are near the 
peak of efficiency, the host halo of a cluster will always be 
significantly above the peak. Thus, this scale cannot probe 
the mismatch at peak efficiency. As larger lensing samples at 
lower masses become available, however, it may be possible 
to perform direct comparisons at those scales. 

The disagreement between abundance matching and profile 
matching could be further influenced by the use of a lumi- 
nosity function instead of the A/, mass function. Because 
the luminosity function varies depending both on galaxy type 
(and thus, color) and choice of band, it could bias the inferred 
abundance matching scales differentially for different galaxy 
types. This explanation for the difference in Figure[TO]is sup- 
ported by results such as lMoster et al.l (l2009h that find a char- 
acteristic scale in the A/vir — A/* relation at A/vir — lO^^A/©, 
just where our profile matching efficiency is highest. 

Other issues affect our interpretation of the dSph galaxies 
in our sample. First, almost all of them are located within the 
virial radius of the Milky Way, meaning that their dark matter 
halos are subhalos, which may follow different scaling rela- 
tions. We consider the effect of this on our derived relations 
in the next section. Also, for the very faintest galaxies, we are 
approaching a regime where surface brightness effects could 
lead to an observational bias to detect only the highest A/vir //> 
galaxies (.Bullock et al., 2 lOl) . 
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Fig. 10. — Comparison between profile matching and abundance matching. The plots show averaged relations for virial mass-to-light-ratio against virial mass 
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(see il5.2t . The gray shaded region is the expected mass-to-light-ratio for a system in which all of its baryons are converted into stars, computed as Tt:/fha.i-yon- 
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Despite these caveats. Figure [TO] does clearly show similar 
patterns to those noted in Sj4] On the faint end, Tvir = M^ij-fL 
shows a much steeper dependence on dark matter mass (this 
time A/vir), while the CSph on the bright end are much more 
sensitive to L. This continues to suggest the dark matter ha- 
los are of greater importance for dSphs, while Es and CSph 
scalings are more controlled by baryonic physics. 

A final intriguing property of the profile matching scheme 
is that there is a built-in consistency check for monotonic- 
ity in the Afvir — L relation. Specifically, if the A/1/2 vs. 
ri/2 relation is anywhere shallower than the Afhaio(?') profile 
it is matching, then the assumption of a monotonic, one-to- 
one mapping from averaged halo mass (and density profile) 
to averaged galaxy luminosity must break down. The fact that 
the model used here does not have this problem implies self- 
consistency, although it does not guarantee this property in 
the actual universe. Clearly, given the size of the measure- 
ment errors (see below) the data at this point are not accu- 
rate enough to determine whether or not the relation becomes 
shallow enough to make the mapping double valued over a 
small r range. We note, however, that if we only consider the 
smallest (magenta, dSph) galaxies (ri/2 ^ 1 kpc), the relation 
appears consistent with Afi/2 oc rf,^. For r << rg (true for 
most of the dSphs here), NFW halos obey A/haio oc r'^12, so 
the profile matching is just at t he limit of monotonicity in the 
relevant halo mass range (see, IWalker et al.ll2009t IWolf et all 
|2pi0, for related discussions). We return to this issue in the 
next section. 

5.2. Uncertainty and Scatter in the Afvir ~ L relation. 

Profile matching to the fundamental curve provides a po- 
tentially strong constraint on galaxy formation models, and 



in principle this method provides a means to test whether or 
not there is an average, monotonic L — A/vir relation between 
galaxy luminosity an halo mass, and to investigate the degree 
of scatter about this relation. Unfortunately, this level of pre- 
cision testing is h indered by several uncertainties. First, as 
discussed in SJ43] there is observational uncertainty that af- 
fects our ability to measure the scatter about and underlying 
shape of the fundamental curve. Second, there is theoretical 
uncertainty in the average mapping between an inner mass 
^'-^u^ ~ -^haio(?'i/2) ^nd halo virial mass, which is particu- 
larly difficult (and somewhat ill defined) for the dSph popula- 
tion we consider because they are subhalos. Finally, even in 
the limit where the theoretical mapping between the average 
M{r) profile and Afvir is perfect, there is a well-known scatter 
in halo profiles at fixed rnass dj ing 2000; Bullock et al. l2001t 
IWechsler et al.l2002l;lBovlan-Kolchin et al..2009.) and this im- 
poses a limiting cosmic scatter in the map between Afi/2 and 
Afvir- We discuss how all of these issues affect the Afvir — L 
relation in what follows. 

Figures [TT] and [12] provide visual presentations of the ob- 
servational and theoretical uncertainties in the profile match- 
ing relations for Afvir vs. L (left) and the equivalent im- 
plied relations of Mvii/L vs Afvir (middle) and My^^/L vs. 
L (right). Starting with observational uncertainties, the error 
bar on Afvir for each data point is estimated by offsetting the 
observables by their la errors in Af-R^ and ri/2, and per- 
forming the profile matching for each data point individually 
(using the mean fundamental curve relation for Model 3). For 
the Es, we use the error bars adopted in the previous section 
(factor of 2 on M^,). We note that this this implies very large 
errors on AfR^ for the E (and dE) galaxies, because these are 
the systems for which Af,/2 is closest to Afi/2, and hence 
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the possible error in M, has the largest effect on M^J^- This 
large uncertainty in MR^ maps to an even larger (relative) 

uncertainty in 7\/vir- Figures [TTI and [T2] are distinguished by 
use of the dMRL-2 and dMRL-3 models, respectively. 

Next we consider the cosmological scatter in the dark mat- 
ter mass enclosed within a given radius for an ensemble of 
halos with identical virial masses. For field halos, this scat- 
ter can be accounted for by the scatter in the concentration- 
mass relation for dark halos, which is approximately log- 
normal in concentrati on with a variance of Alog(c) = 0.14 
dWechsler et al.ll2002 !). In principle, this cosmic scatter pro- 
vides a lower limit on point-to-point scatter that can be mea- 
sured in a profile matched A/vir — L relation. We illustrate the 
magnitude of this cosmic scatter by the middle (dark gray) 
shaded band, which traces our best-fit relation (shown as a 
solid blue line connecting blue circles) in each panel. We 
see that this cosmic variance is particularly important for the 
smallest galaxies. This cosmic variance scatter is the mini- 
mal possible scatter expected for galaxies in ACDM. Even if 
galaxy properties tracked virial mass in a precisely one-to-one 
fashion, they would scatter about the profile matching relation 
with at least this amplitude. Q 

An additional component of scatter and uncertainty must 
be considered for the dSph galaxies - because they are satel- 
lites of the MW, their dark matter halos are subhalos, and 
hen ce do not obey the s ame scaling re lations as field halos 
(e.g lBuUock eTani200!t ISpringel et al.ll2008ft . More specifi- 
cally, it is inappropriate to speak of a virial mass for a subhalo, 
because subhalos tend to be tidally truncated at radii that are 
smaller than the virial radius they had when they were first 
accreted. A more meaningful mass to be associated with each 
dSph is its halo's virial mass at the time it was accreted. It is 
this mass, the virial mass at accretion, that would most likely 
show a strong correlation with galaxy luminosity. 

Two competing effects may act to modify standard (field) 
mapping between inner mass and virial mass. First, at fixed 
virial mass, a halo at higher redshift will tend to be denser 
at a fixed physical radius than a halo of the same virial mass 
at a later redshift (because the virial density scales roughly 
with the density of the universe). Therefore, if a subhalo was 
accreted at some high redshift (e.g. z — i) and it experienced 
no mass loss in its central regions (unlikely) then our virial 
mass estimates are biased high. The lower (red) shaded region 
in the L < 10^ L© band of Figure[TT|illustrates the degree by 
which the median relation would need to be shifted down in 
order to account for a z < 3 accretion that experienced no 
mass loss within its central region after accretion. The lower 
edge of the red band corresponds to the relation expected if 
all dSphs were accreted at z = 3 with no mass loss. 

The second, competing processes that adds uncertainty 
to profile matching estimates for subhalos is tidal mass 
loss. Halos tend to lose mass at all radii after they are 
accreted, and this acts to decrease their central density for 
a fixed vir ial mass at accretio n . The cosmological simu- 
lation of Bovlan-Kolchin et al.l (I2009h shows that the me- 
dian subhalo at z = in a Milky- Way-type host has lost 
75% of its initial total mass, while ~ 85% of subhalos 
in have lost < 90% of their initial total mass (Boylan- 
Kolchin 2010, private communication). However, the mass 



' In principle, if galaxy luminosity had a secondary dependence on halo 
concentration, then the covariance could act to reduce the cosmic scatter from 
profile matching, but this seems tuned and unlikely. 



loss is far less signi ficant in the inner regions we are prob- 
ing here {K azantzid is et all 120041: iPenarrubia et al.l l2008cl: 
IWetzel & White 2010; 'PenarrubiaeialJ|2010|). The simula- 
tions of Bullock & Johnston (2005) show that a 75% (90%) 
loss of total mass, results in a mass loss fraction within the 
inner 300 pc of only 20% (40%) - where ri/2 = 300 pc is 
the median half-light radius for our d Sph sample. For th e 
mass range of relevance Afvir oc M^q^ (iBuUock et al.ll2010l) . 
which implies that our fiducial Afvir determination from field 
halo profile matching would be under-estimated by a factor of 
(0.8)"'^'^ ^ 2 for median subhalo mass loss, and by a factor 
of (0.6)^'^'^ ^ 5 in the case of 90% total mass loss. Thus, 
in Figure [TT] we include an upper (green) shaded region cor- 
responding to a factor of 5 increase in the inferred A/vir, as a 
conservative estimate of the maximal scatter. This treatment 
is conservative because we expect that systems with the most 
mass loss will also have been accreted earlier, and therefore to 
have had higher virial densities overall. This offsetting effect 
has been ignored in the upper green shaded band. 

Of course, if we knew the redshift of accretion and orbital 
trajectory (including mass loss) of each dSph in our sample, 
we could perform the profile matching in a more exacting 
way, but this is not practical with present d ata due to the un- 
certainties in the orbits of the MW satellites (iLux et al. 2010). 
Therefore, we have added both the effects in quadrature to the 
cosmic variance error band in Figure [TT] in order to derive a 
limiting scatter estimate shown as the outer light gray regions 
in Figure [TTI 

Thus, the shaded bands about the average relations in Fig- 
ure[TT]can be thought of as a limiting theoretical scatter about 
the relation. In principle, if the data at a particular scale scatter 
about the relation with a larger variance than indicated by the 
shaded band, then this would be indicative of intrinsic scatter 
in the Afvir — L relationship. This then implies that the sec- 
ondary scalings in galaxy formation (e.g. 2-D scalings such 
as the fundamental plane) can be fit to the data set to pro- 
vide useful information. Conversely, if the scatter (including 
observational errors) is consistent with the theoretical scat- 
ter, the secondary scaling relations cannot be measured at that 
scale. 

The possibility of detectable intrinsic scatter is is partic- 
ularly interesting at the faint end, where it has been noted 
that despite the wide ranges of luminosities, the MW dSphs 
appear to have similar halo masses, albeit with large scatter 
(Strigari et al. 2008a; Wolf et al. 2010) - this could be due to 
scatter (observational or intrinsic) masking a weak relation, 
scatter in halo mass about a new scale in galaxy formation, 
selection effects (e.g. the stealth galaxies' influence as dis- 
cussed in S]5j, or some as yet unknown effect. This scale ap- 
pears particularly strongly in[T2]due to the preferential fitting 
on the dSph. These data also admit a steepening power law 
instead of a true scale at the low-mass regime as suggested by 
iKravtsovl (|2010) to match the dSph luminosity function, so 
we plot this relation in Figure [TOl 

We note in Figure [TT] that there is a systematic offset for 
the bright dSphs. This is primarily due to the tension between 
fitting the RM relation for the dE and the bright dSphs with 
a single power law, as is used for dMRL-2. In Figure [12] 
this offset is essentially gone, as the fit in the RM relation is 
tailored to fit best for the dSphs. This comes at the price of a 
poorer fit for the other galaxies, however, as well as an anoma- 
lously low Afvir/ i apparent in Figure [TO] (green-dashed line 
in lower-right panel). It is unclear if this tension is due to 
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Fig. 1 1. — Point-to-point scatter about tlie median A/vir — L relations: V-band luminosity against virial mass (left), virial mass-to-light ratio against virial 
mass (middle), and virial mass-to-light ratio against luminosity (right). The blue line with points is the median relation for the dMRL-2 model. The width of the 
shaded band is the minimum galaxy-to-galaxy scatter that we expect in the limit where virial mass con'elates with L with zero scatter. The band is not the en'or 
in the A/vir ~ L relation, but rather the minimum derived scatter set by a fundamental cosmological uncertainty in the mapping between M{r--ii2) = ^^^12 ^^'^ 
Afvir. On the faint end, where the galaxies in our sample are in sub-halos, the inferred mass can be scattered upwards due to uncertainty in infall redshift (green 
inner shaded region), and downward due to mass loss (red inner shaded region). Adding these effects in quadrature with the cosmic variance (middle dark gray 
shaded region) results in the full cosmological scatter (outer Ught gray shaded region). Points represent the objects described in ijl] with virial masses inferred 
by inverting the grid of Figure[8]and error bars due to the effect of the observational error bars on the mapping, including a factor of 2 uncertainty in A/* for dEs 
and Es. For clarity, eiTor bars are only shown for the 20 E galaxy bins with the largest number of galaxies. The MPJ"^ assumed for the Es has been connected 

for baryonic contraction using the estimate described in ij5] For the CSphs, it is important to note that the vertical error bars are likely underestimated due to 
the uncertain effect of using the cluster galaxies instead of the intra-halo Ught in measuring the velocity dispersion. The data point color and point-type scheme 
matches that of the previous figures (see JS). Intrinsic scatter in the Afvir ~ L relation can only be detected in cases where the observational error bars are small 
enough to detect a dispersion smaller than the shaded band. We see that the population of dSph galaxies is close to this level now. 
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Fig. 12. — Same as Figure [771 but using the dMRL-3 model to compute the profile matched relation (blue solid line). The nature of the "bump" feature is 
discussed in Js] 



problems with ACDM accounting for the existence of galax- 
ies in the halos of the bright dSph, evolutionary effects on 
subhalos (as discussed above), or the influence of baryonic 
contamination of Mi/2> which is unaccounted for in our anal- 
ysis of the dSphs. 

Unfortunately, as is clear from comparing the data points to 
the shaded band in Figure[TTl the observational uncertainty is 
still slightly too large on the faint end to determine if there is 
significant intrinsic scatter about the fundamental curve, al- 
though the data are close. Similarly, while the very faintest 
galaxies show deviation from fundamental curve in a way 
consistent with a new scale of flat M^i^jL, this level of devia- 
tion is not statistically significant. Similar uncertainties likely 
apply to M3 1 satellites, making it difficult to interpret the pos- 
sible existence of an offset (Kalirai et al. 2010) . Marginal 
improvements in data quality may be enough to shed light 



on these questions, however, as observational errors could be 
brought to the level of cosmological scatter Furthermore, 
the predicted existence of far more faint dS phs in the Lo- 
cal Group to be detected in upcoming surveys dToUerud et alj 
l2008t iMartin et al.ll2009l: iBullock et alJIIoTol) provides hope 
that this degeneracy between intrinsic and observational scat- 
ter may be broken by sheer statistics. Nevertheless, the cur- 
rent data are not good enough to definitively address this ques- 
tion. 

There is also hope on the bright end. Interestingly, the most 
massive, luminous objects are the ones that face the least cos- 
mological scatter associated with the profile matching tech- 
nique. As can be seen in Figure[8] as one travels along the fun- 
damental curve projection to large values of ri/2 ^nd MR^, 
the associated il/vii determinations become more cleanly de- 
fined. Unfortunately, it is in this regime where our inabil- 



Fundamental Curve 



19 



ity to determine CSph velocity dispersions limit the ability to 
cleanly determine Mi/2- 

6. CONCLUSIONS 

We have examined the scaling relations for a broad collec- 
tion of spheroidal stellar systems in an intrinsic MRL space of 
half-light mass (Afi/2; Equation [U, half-light radius (ri/2), 
and half-luminosity (Li/2). These MRL coordinates are a 
theoretically-motivated transformation of the familiar funda- 
mental plane variables and can serve as a bridge between di- 
rect observables and the predictions of galaxy formation mod- 
els. The latter is facilitated by considering an alternative space 
we refer to as dMRL space. In dMRL space, the mass vari- 
able is M^J^ - the dark matter mass within rij2 rather than 

the dynamical mass. 
Our main findings are as follows. 

1 . All spheroidal galaxies — stellar systems with their own 
dark matter halos — track a 1-D fundamental curve 
through MRL space. This curve is visualized in 3-D 
in Figure |6] and represented analytically in Equations 
|4] and |5] (with parameters from Table |2|i. The fun- 
damental mass-radius-luminosity relation transitions 
from Ml/2 cxi r^f2 oc L'^f2 for the faintest dwarf 

spheroidal (dSph) galaxies to A/1/2 '^ ''1/2^ '^ ^ifi 
for the most luminous cluster spheroids (CSphs). This 
ri/2 — Li/2 scaling (MRL-2 model) is a good fit for the 
dSphs if we take into account the fact that the lowest 
luminosity dwarf galaxies suffer from surface bright- 
ness incompleteness (which biases the sample towards 
smaller ri/2)- If we ignore this bias, the raw empirical 
relation (MRL-1 model) gives A/1/2 oc r\f^ oc L^^?'^ 
on the faint end. 

2. Dwarf ellipticals (dEs) and normal ellipticals (Es) in- 
habit the transition regime between the limiting power 
laws, where the dynamical mass-to-light ratio within 
ri/2 is minimized at T1/2 ~ 3. The dynamical 
mass as a function of ri/2 transitions quite abruptly 
as the galaxies become baryon-dominated (see Figure 
|7J. When we subtract out the baryonic component 
with estimates for the stellar mass (although these are 
subject to uncertain systematic errors), the relation is 
better fit by a power law, particularly when we in- 
clude an estimate for the effect of baryonic contraction 
(see the inset of Figure [8]l. The inferred slope for the 
ri/2 - A/R^ relation is A/^^ oc r\f^, sUghtly steeper 

than the M (xr^ relation that has been discussed in the 
literature (Gentile et all I2009t iN^olitano et al] 120 lOt 

IWalkeretal. 2010). 

3. Globular clusters (GCs) and ultra-compact dwarfs 
(UCDs) do not follow the fundamental curve relation. 
Instead, GCs and UCDs inhabit overlapping/connecting 
regions in MRL space that resemble sections of mass- 
follows-light planes near A/1/2 = 3 /^i/2, as illustrated 
in Figure m See Equation |2] for the exact form of the 
plane that separates this GC locus from the dSph por- 
tion of the fundamental curve. Note that the UCDs in 
our sample exhibit a small "tilt" away from the mass- 
follows-light plane, while GCs exhibit no such tilt - 
thus it cannot be ruled out that UCDs are a part of 



the galaxy sequence, but are intrinsically rare in the re- 
gion where they meet the fundamental curve. However, 
dSphs separate distinctly from GCs and UCDs in MRL 
space, implying that if UCDs are actually embedded in 
dark matter halos, an irreducible dichotomy exists in 
galaxy formation. 

4. The fundamental curve relation in dMRL space allows 
us to connect galaxies to their dark matter halos via an 
approach we call profile matching. Specifically, at each 
luminosity, an average galaxy sits in a specific point in 
the A/P.^ — ri/2 plane. This mass-density point can 
be mapped to an average dark matter halo virial mass, 
as illustrated in Figure [8] While we assume standard 
NEW halos in ACDM, this technique is easily adapt- 
able to any dark matter halo type that can be cast as a 
single-parameter family. In the end, we can construct 
relationships between luminous galaxy properties and 
their dark matter halo masses. This profile matching 
technique for deriving the A/vir — L is most accurate at 
the high and low luminosity extremes (where dark mat- 
ter fractions are highest) and is therefore quite comple- 
mentary to statistical approaches that rely on having a 
well-sampled luminosity function. 

5. Independent of any global abundance or clustering in- 
formation, we find that (spheroidal) galaxy formation 
needs to be most efficient in ACDM halos of virial 
mass A/vir — 10^^ A/0 and to become sharply ineffi- 
cient in masses smaller than A/vh- < 10^*^ Mq. On the 
other hand, the inefficiency of galaxy formation seems 
occur more gradually as halos become more massive 
than A/vir — 10^"^ A/0. Rather, the inefficiency sets 
in sharply in luminosity at L ~ 10^^ L0. These re- 
sults are qualitatively consistent with the expectations 
of abundance matching (see Figure [TO] i, although only 
if we use models that account for surface brightness se- 
lection effects on the faint end (dMRL-2 and dMRL-3). 
The sharpness of the transition with A/vir on the faint 
end may imply the dark matter halo or potential depth 
drives scaling relations for low-mass galaxies, while the 
stronger dependence on L on the bright end suggests 
baryonic physics controls the massive galaxy regime. 

6. Object-by-object scatter about the A/vir — L relation re- 
mains very difficult to quantify. Nevertheless, despite 
the large theoretical uncertainties associated with our 
profile matching technique at the low-mass end, the ob- 
servational data for dSphs are almost to the point where 
we can explore intrinsic scatter about this relation in 
the smallest systems. On the other hand, the theoret- 
ical uncertainty in the mapping between points in the 
M^J^ ~ '^1/2 plane and halo virial is much smaller on 
the scale of CSphs, so there is promise at the bright 
end from that respect. Unfortunately, stellar velocity 
dispersion for CSphs remain very difficult to obtain di- 
rectly. A better approach would be to consider alterna- 
tive mass-radius measurements for CSphs (based, for 
example, on X-ray studies) as has r ecently been ex- 
plored bv lTruiillo-Gomez et alj (I2010I) . 

We close by mentioning that the existence of a fundamental 
curve in MRL space is not out of line with an understanding 
that galaxy properties show strong correlation with a single 
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parameter (see, e.g. iDisnev et alj|2008l for similar results on 
an Hl-selected sample). Nevertheless, this fact does not imply 
that all galaxies belonging to a given evolutionary sequence 
are completely or even primarily controlled by a single pa- 
rameter - only that their first-order scaling relation is charac- 
terized by a single parameter when the galaxy properties are 
averaged. Our viewpoint is rather that the MRL relation pre- 
sented above provides a useful bridge between observational 
properties and theoretical models. At the very least, mod- 
els should be able to reproduce the 1-D scaling relation pre- 



sented. Some guidance to that aim is provided by our dMRL- 
inspired profile matching, which seeks to unite galaxies across 
a space of virial mass, stellar luminosity, and stellar radius 
self-consistently. 
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APPENDIX 
ANALYTIC FIT TO HALO MATCHING RELATIONS 



In order to provide an analytic description of the derived My-„ relations for our fiducial (dMRL-2) model results, we perform a 
least-squares fit of y vs. Afvir for each of y = Li/27 ^1/2^ Afi/2' ^nd (Mvii/L) using the same fitting form as Equation|4] 



log 



M± 



B 



S-MiA-B) 



arctan(A^/M^) 



(Al) 



Here, Ai = log(Afvir/Af°ii.) defines a characteristic virial mass scale M°jj. at y = j/o and W sets the width of the transition from 
y ex M^j. and y oc M^j. at small and large Afvir, respectively, and S sets the offset in log(j/) over the transition region. Each of 
the fit parameters A, B, W, S, M^^^., and j/o ^e provided in Table |3]for our four y relations, corresponding to the four panels of 
Figure |9l 

We find that the A/R^ vs. A/vii and ri/2 vs. Afvir relations for dMRL-2 are also very well characterized by power-laws. 
Specifically we find 



and 
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The i-to- Afvir relation, meanwhile, can be approximated on the faint end as L ex Af.^;-^^^ and flattens on the bright end to 



L oc Af Oj^e 



ALTERNATIVE DATA PROJECTIONS 



The fundamental plane of bright elliptical galaxies (iDiorgovski & Davislll987l : lDresslerlll987l : iFaber et alj|1987b lies within 
a 3-D parameter space that consists of the velocity dispersion (a), the 2-D half-light (effective) radius (Re), and the surface 
brightness (fg). Here we define fg such that it is the mean surface brightness within f?e, in units of Lq pc~^, although we note 
that slightly different definitions are sometimes used in the literature. The fact that these three variables are direct observables 
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TABLE 3 
Profile Matching Model Parameters for Eouation IAII 
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that scale together motivates the consideration of galaxies in this space. Two-dimensional projections of our data set (Table [U 
on the fundamental plane axes are shown in the three panels of Figure[T3] In Figure[14]we plot in a 3-D rendering of these same 
data. Also shown in transparent red in Figure[T4]is the best-fit fundamental plane of Graves et al. ( 2009g). 

From these plots, it is apparent that while the normal elliptical galaxies (E) lie well within the fundamental plane of l&aves et al] 
(l2009b ). the CSphs and dS phs lift away from th e plane in a non-trivial manner, in contrast to data sets that do not reach those 
extremes in luminosity (e.g. lBurstein et al.ll997l) . How ever, it has been noted in the literature that the f aint end of the fundamental 
plane (towards dEs) shows curvature up off the plane jZaritsky e t al."20 06al: iHvde & Bernardill2009h . and bright-end deviations 
from the fundamental plane are discussed in further detail inlZaritsky e t al.l ( I2006bl) . Here we note that the separation from 
the plane is much more significant when the dSph galaxies discovered in the SDSS (Re < 450 pc) are included alongside the 
"classical" dwarfs, as the SDSS dSphs extend nearly perpendicularly from the fundamental plane, with that in mind, the deviation 
from the plane is significant far beyond the scatter in the fundamental plane derived for bright E galaxies. 

The "tilt" of the E fundamental plane here can be interpreted in the context of this curvature; the tilt of the fundamental plane 
is simply the shift of the observational fundamental plane from the expected virial plane (see Sj3] and the two planes in Figure 
[3]). Curvature off the plane is then just continuation of this tilt past the typical regime of Es. The tilt in the Es can potentially 
be driven by a mix of stellar m ass-to-light ratio variations and/or vari ation in the dark raatter-to-baryon fraction within the halo 
of the galaxy in q uestion (Cap pellari et alJl2006l:lBolton et a l. 2007; Humphrey & Buotell2010l;lNapolitano et alJl2010l;lTreu et al] 
I20 lot [Graves & Faber.2010 ). It may also be an aperture affect due to dissipation caus ing a chang e in the app arent dark matter 
fraction by packing more baryonic material in the same volume of dark matter halo (Ro bertson et al. 2006; Hop kins et alJ lTOOS). 
For our purposes, however, it is sufficient to note that the magnitude by which the largest and smallest spheroidal galaxies peel 
away from the fundamental plane cannot be explained by baryonic effects - it can only be addressed in terms of dark matter 
content, due to the very large mass-to-light-ratios. 

For comparison with other work. Figure [15] shows the projection of the data onto the kappa (k) space of iBender et al.l (Il992h . 
a coordinate rotation that enables a reasonably physical interpretation with ki oc logM, K2 oc log [{M/L)lf), and K3 oc 
log(M/L) such that ki and K2 define a plane that is approximately parallel to the fundamental plane for ellipticals. 
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Fig. 13. — Two dimensional projections of the three-dimensional data shown in Figure [T4l along the coordinate axes. The color/symbol code maps to galaxy 
type as indicated, matching the scheme of Figure[T] 



In Fi gure [T6l we show this data set again in MRL space (as in Figure O, but we no w overplot the fundamental manifold of 
IZaritskv et al. (2006bj) in transparent green. We use the fundamental manifold from Zaritskv et al.l ( 12008.) Table 1, using the 
transforms from jS]t o convert from fundamental plane space to MRL space. We also adjust the luminosity fro m the value for 
I-band (used in iZaritskv et al.112008 ) to V-band assuming all objects have V-I colors of typical E galaxies from lFukugita et aT] 
(11995). While there will be an additional bias because Re for V and I bands will differ, this is likely small relative to the scatter 
and hence we disregard it. 

Finally, we show the mean mass density of our data set as derived from the middle panel of Figure |2] This derived simply by 
tilting the log(ri/2) — log(A/i/2) relation to give density within ri/2 instead of mass (i.e. residuals from the dashed-dotted line 
in the middle panel of Figure|2l but with a different normalization fixed to standard density units). 
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Fig. 14. — Three dimensional representation of the data set in fundamental plane coordinates of log(/e/i0Pc~^), log(i?e/kpcV and log((T/km/s). The 
red transparent plane is the fundamental plane for elHpticals from Graves et al. (2009b), and the color/symbol code matches Figure[T] A rotating animation of 
this plot is available at,http : //www, physics . uci . edu/ ~bul lock /f curve /movies . html, 
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Fig. 15. — Two dimensional projection of the three-dimensional data shown in Figure [T4l onto the k space of IBender et al.l jl992f) : ki oc log A/, K2 oc 
log [{M/L)!^), and Ka oc log (M/ L), and ki and K2 are approximately parallel to the fundamental plane. The color/symbol code maps to galaxy type as 
indicated, matching the scheme of Figure[T] 
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Fig. 16. — Three dimensional representation of the fundamental manifold of lZaritsk vetalJ J2006bl) i n MRL space (i.e. Figure |3). The transparent (green) 
manifold with solid lines is the fundamental manifold with the fit coefficients from lZaritskv et Jl <20Q8D corrected to V-band (see text). The (blue) plane with 
dashed borders corresponds to Afj /2 = SL^ /2, the mass-follows-light plane. The data point color and point-type scheme matches that of the Figure[3](or see 
^. A rotating animation of this plot is available at.http: //www .physics .uci ■ edu/ ~bul lock /f curve /mo vies .htmll 
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Fig. 17. — Mean mass density within rj /2 as derived from the r^ /2 — M^ /2 relation (e.g. middle panel of Figure|2). The left panel uses the raw dynamical 
mass A'Ii/2 to compute the density, while the right panel uses the mass AfPi* for which the stellar contribution has been subtracted. To guide the eye, we include 

the Pi/2 oc r~^ relation (black dashed lines), and for the right panel, the isothermal case of Pi/2 oc t~'^ (black dotted line). The color/symbol code maps to 
galaxy type as indicated, matching the scheme of Figure [T] 
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UCD MASS ESTIMATES 



As discussed in ij3] UCDs present a puzzle in the MRL space. The most massive UCDs approach the fundamental curve, 
although with a gap that could potentially be a result of selection effects. Regardless, the current sample of UCDs form a distinct 
group (with GCs) from the dSphs for the faintest/smallest objects. Thus, if this sample of UCDs and dSphs are both galaxies, the 
scaling relations split into a dichotomy or bimodality at the faint end. For this paper we have focused on the dSph side of this 
relation, but here we consider the UCDs in the profile matching context. 

Figure [18] is analogous to Figure |8] but zoomed in on the faint end and with UCDs added. Note that for UCDs, we determine 
the dark matter mass for the UCDs as M^J^ = M1/2 — ^1/2"^, where an T is taken to be fixed at 1 (open circles) or 2 (filled 
circles). For the latter, we also show error bars based on a possible factor of 2 systematic uncertainty in stellar models, based on 
the discussion in il4.3l for E galaxies. The dSph error bars are from|8] based on the observational error bars in ri/2 and a. 

The grid of NFW halos in Figure [18] clearly shows that the implied dark matter densities for UCDs are most consistent with 
cluster-sized (or larger) dark matter halos (Afvir ^ 10^^) . Taken at face value, this is impossible, as there are not enough of such 
halos where UCDs are found, and they would have clear kinematic effects on neighbors if UCDs had such massive halos. A few 
possibilities might explain these large virial masses. If UCDs do indeed have dark matter halos, baryonic contraction might boost 
their central densities (as described in i2]for E galaxies). However, given the extreme stellar densities and small sizes of UCDs 
(and hence short dynamical times), it seems unlikely that any baryonic contraction would be adiabatic. Hence we cannot apply 
adiabatic contraction as we have used to correct masses for the Es. A baryonic contraction model appropriate for UCDs could be 
used in the same way, although we do not do such a correction as such a model does not yet exist. 

An alternative possibility is that the stellar population estimates are systematically in error. The error bars in[T8]imply that such 
errors could explain most (possibly all) UCDs as entirely stellar objects-this corresponds to those where the error bars are upper 
limits. Alternatively, they may have dark matter halos with much smaller virial masses, but without better stellar population 
models, there is no way to tell the difference. 

Thus, while UCDs are possibly consistent with lying inside dark matter halos, this implies a dichotomy in galaxy formation as 
well as being impossible to explain with standard ACDM dark matter halos. We therefore favor the simplest view that they are 
purely stellar systems. 
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Fig. 18. — Profile matching schematic for UCDs and dSphs, based on Figure[8] Open blue circles are UCDs with mass-to-light ratios (T) of 1 assumed for 
computing AfRj . while solid blue circles assume T = 2, and error mass error bars assume a factor of 2 uncertainty in T. Upper limit error bars (extending 
to the bottom of the plot) correspond to those where the systematic uncertainty could result in a UCD fully consistent with its stellar mass implied from the 
luminosity. Magenta triangles are dSph galaxies, with observational error bars. The colored lines are the grid of Figure [8] representing the mass profiles for a 
range of NFW halos. Thus the densities of UCDs imply that they are either within the most-massive halos, or primarily purely stellar. 



